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Abstract: Within gauge/gravity duality, we consider finite density systems in a heli¬ 
cal lattice dual to asymptotically anti-de Sitter space-times with Bianchi VII symmetry. 
These systems can become an anisotropic insulator in one direction while retaining metallic 
behavior in others. To this model, we add a 17(1) charged scalar and show that below a 
critical temperature, it forms a spatially homogeneous condensate that restores isotropy 
in a new superconducting ground state. We determine the phase diagram in terms of the 
helix parameters and perform a stability analysis on its IR hxed point corresponding to a 
hnite density condensed phase at zero temperature. Moreover, by analyzing fluctuations 
about the gravity background, we study the optical conductivity. Due to the lattice, this 
model provides an example for a holographic insulator-superfluid transition in which there 
is no unrealistic delta-function peak in the normal phase DC conductivity. Our results 
suggest that in the zero temperature limit, all degrees of freedom present in the normal 
phase condense. This, together with the breaking of translation invariance, has implica¬ 
tions for Homes’ and Uemuras’s relations. This is of relevance for applications to real world 
condensed matter systems. We hnd a range of parameters in this system where Homes’ 
relation holds. 
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1 Introduction 


Significant progress has recently been achieved in applying gauge/gravity duality to strongly 
coupled systems of relevance to condensed matter physics. In particular, different ap¬ 
proaches were proposed to include a lattice into the dual gravity background, in order to 
holographically study the conductivity in systems with broken translation invariance. Sys¬ 
tems with manifest translation invariance display an unrealistic h-function at zero-frequency 
in the conductivity; breaking the symmetry weakly broadens this into a realistic Drude peak 
known from condensed matter physics as a consequence of momentum dissipation. Hologra¬ 
phy allows moreover the exploration of the consequences of translational symmetry breaking 
for strongly correlated systems beyond the Drude peak, both in the weak Drude regime and 
for stronger lattice potentials. We will follow this avenue in the present paper. 

Within holography, translation breaking approaches include explicit breaking by a mod¬ 
ulated Ansatz for the chemical potential [1-4], the use of massive gravity [5, 6] or linear 
axions [7-9], or other lattice Ansatze such as the Q-lattices [10] or the method used in this 
work, Bianchi symmetric solutions [11, 12]. In some cases, translation invariance is also 
spontaneously broken, for instance when Chern-Simons oi F A F terms are present in the 
gravity action [13-20], or in an external SU{2) magnetic field [21, 22]. Explicit breaking 
with an interesting phenomenological consequence is realized in the helical lattice approach 
[11, 12, 23-27]. The original motivation to study this model was that the helical symmetry 
allows for momentum relaxation along one spatial direction without the need to solve PDEs. 
The helix in one of the field theory directions along the boundary is encoded in a non-trivial 
background 17(1) gauge field on the gravity side of the holographic duality, and its shape 
is protected by a so-called Bianchi VIIq symmetry. In addition to that ‘helix 17(1)’, the 
five-dimensional gravity action (2.1) that we study involves a separate ‘charge 17(1)’ dual to 
a globally conserved charge current in the boundary theory. This is needed to encode a field 
theory at finite density. As discussed in [26], the natural finite density state of this model 
is a conducting metal, but it can display a transition to an insulating phase as a function 
of the helix momentum. The remarkable aspect is that this new phase is uni-directional 
and anisotropic: It is an insulator only along the direction of broken translation invariance 
along the helical axis. In the orthogonal directions the system remains a metal. From a 
condensed matter point of view, this system resembles a so-called quantum smectic. 

The specific objective we shall be interested in this paper is the consequences of transla¬ 
tional symmetry breaking for the transition to superconductivity. This was also recently 
studied in a holographic Q-lattice in [10, 28] and in axion and related holographic super¬ 
conductor models in [29, 30]. In both cases only isotropic models were considered, though 
both models can support anisotropic lattices [31, 32]. In our intrinsically anisotropic helical 
lattice, the dual gravitational dynamics imply that the favored ground state will neverthe¬ 
less be an isotropic s-wave superconductor. In fact the only Bianchi VIIq symmetric and 
time-independent Ansatz for the scalar field dual to the order parameter is just a constant 
in boundary direction. This is the system we shall study. We show that a scalar field added 
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to the helical lattice action and charged under the second U{1) gauge field condenses below 
a critical temperature, both in the insulating and in the conducting phase. We explore the 
phase diagram which is determined by the amplitude and the momentum of the transla- 
tionally symmetry breaking helix, both at finite and at vanishing temperature. Moreover, 
we analyse the thermodynamical as well as transport properties of the different phases — 
metallic, insulating, condensed — of the helical lattice model, and obtain the finite tem¬ 
perature phase diagram of the system. 

Our findings can be summarized as follows: 

The superconducting phase transition: In Section 2 we investigate the finite temperature 
phase diagram of our model, which displays a second order mean field transition from both 
the insulating as well as metallic phase to a superfluid phase. We in particular show that 
the critical temperature Tc depends strongly on the amplitude of the helix, but to a first 
approximation rather weakly on its momentum. This indicates that the strength of the 
translational symmetry breaking (the depth of the potential wells) is more important than 
their spatial distribution. For large amplitude, Tc can in principle be suppressed all the 
way to zero, and a quantum phase transition to the uncondensed phase may be expected. 

The critical temperature Tc does depend mildly on the helix momentum p, in a curious 
way. Starting in the phase at small helix momentum which is originally a zero-temperature 
insulator in the normal phase, Tc decreases with increasing helix momentum. However, Tc 
grows again for even larger helix momentum. This might be understood from the obser¬ 
vation that initially with increasing momentum the underlying original insulating system 
changes to a zero temperature conductor in the normal phase, but then for even larger 
momentum turns back into an insulator. 

The optical conductivity: In Section 3 we calculate the optical conductivity in the direc¬ 
tion of translational symmetry breaking in the insulating, conducting as well as condensed 
phases. In the insulating and conducting phases we reproduce the results of [26]. In the 
condensed phase we observe the appearance of a gap at low frequencies as expected for 
spontaneous symmetry breaking. The spectral weight is transferred to a (5-function contri¬ 
bution at zero frequency: this is confirmed with the Ferrell-Glover-Tinkham sum rule. 

The virtue of the helix model is that this (5-function is now cleanly interpreted as the con¬ 
sequence of spontaneous symmetry breaking. There is no artificial contribution due to 
translational symmetry. The strength of this (5-peak therefore defines the superfluid density 
in the condensed phase. For weak momentum relaxation, <C 1 we find that in the limit 
T —7> 0, the superfluid density coincides with the total charge density in the system, as 
measured by the second 17(1) gauge field. This can be understood by the fact that the 
zero-temperature normal state of our system is already in a cohesive phase, in which no 
uncondensed charged degrees of freedom are present in the deep IR. This however does not 
mean that we are dealing with a plain vanilla superconductor. At any finite temperature. 
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the horizon does carry charge. This reflects itself in the temperature dependence of the 
superconducting gap. We And that the low T behavior of the superconducting gap is alge¬ 
braic, i.e. cr{LO*) ~ rather than exponential. Nevertheless as stated earlier, computing 
the thermodynamical charge density Ug and the superfluid density ps independently, we 
And that they coincide in the limit of zero temperature, in the regime of weak momentum 
relaxation. 

The helical system considered has therefore two important properties: Translation symme¬ 
try is broken and all charged degrees of freedom condense at very low temperatures. The 
combination of both these facts enables us to take a further look at Homes’ relation in the 
context of holography. This empirical law, found experimentally [33, 34], states that there 
is a universal behaviour for classes of superconductors that relates the superconducting 
density ps at zero temperature to the DC conductivity (Tdc at 

Ps{T = 0) = Ca-DciTc) ■ Tc- (1.1) 

The constant C, which is dimensionless in suitable units, is experimentally found to be 
around C = 4.4 for in-plane high-Tc superconductors as well as clean BCS superconductors 
and around C = 8.1 for c-axis high-Tc materials and BCS superconductors in the dirty 
limit [33, 34]. Generally, a relation of this type is expected for systems which are Planekian 
dissipators [35]. Homes’ relation was first considered in the context of holography in [36], 
where it was found that for a holographic realization, both translation symmetry breaking 
and the condensation of all charged degrees of freedom are necessary conditions. Both of 
these conditions are realized in the helical lattice system in the present paper. It was found 
in [4] that a simple breaking of translation invariance by a modulated chemical potential is 
not sufficient for a holographic realization of Homes’ relation, essentially since in the limit 
of vanishing chemical potential, the DC conductivity diverges while the superconducting 
density remains finite. Indeed, Homes’ relation cannot hold for weak momentum relaxation. 
However, motivated by the arguments given above, we considered Homes’ relation in the 
context of the helical lattice model for strong momentum relaxation. In a parameter region 
around the minimum of Tc found in Section 2, C appears to be roughly constant for a 
significant region in parameter space. We And a value of about 

6* = 6.2 ±0.3, (1.2) 

which lies between the experimental results for high Tc and dirty limit BCS superconductors. 
These encouraging results call for further detailed analysis, which we leave for future work. 

The zero-temperature ground state: In Section 4 we provide a preliminary analysis of the 
zero-temperature ground states of the condensed system. The starting point is the IR ge¬ 
ometry of the insulating ground state geometry of the helical model in the absence of a 
condensate [26]. We show that this Ansatz for the IR geometry can naturally be extended 
to the superconducting solution. Besides this, the insulating as well as metallic ground 
states of [26] continue to exist. We perform the usual IR fluctuation analysis and delineate 
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IR 


Figure 1. Sketch of the holographic renormalisation group flows without (left panel) and with 
superconducting order parameter (right panel), for the case of large k > 0.57 in which no additional 
unstable bifurcation fixed points appear in the model of [26]. Without the superconducting order 
parameter, the transition from the insulating (red dotted) to the metallic AdS 2 x (blue solid) 
ground states occurs as the helix pitch p/^^ is increased [26]. With the superconducting order 
parameter a new superconducting (green dash-dot-dotted) ground state, (4.1), appears. While the 
metallic ground states becomes dynamically unstable and presumably flows to the superconductor, 
the insulating states stay dynamically stable but are presumably thermodynamically disfavoured 
compared to flows to the superconducting fixed point. 


the various IR relevant directions, if present, and the IR irrelevant directions. This allows 
us to understand the RG flow of the model in principle (Figure 1). We in particular find a 
difference in the instability mechanisms of the insulating and metallic fixed points of [26]: 
While the metallic AdS 2 x becomes dynamically unstable at low temperatures towards 
condensation of the superconducting order parameter, the insulating fixed point stays dy¬ 
namically stable, but most presumably becomes thermodynamically disfavored compared 
with flows to the superconducting fixed point. Since our finite temperature phase diagram 
(Figure 5) indicates the possibility of a zero temperature phase transition between the con¬ 
densed and insulating solution, these results calls for a more detailed study of the phase 
diagram at finite and zero temperature in a future work [37]. 

We conclude in Section 5 where we discuss in particular the implications of our results for 
the holographic realization of Homes’ relation. We end by giving an outlook to further 
investigations. 
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2 Holographic S-Wave Superconductors on a Helical Lattice 

In this section we first explain our setup, which is based on the model of [26]. We then 
discuss and present our results for the finite temperature phase diagram. 


2.1 Holographic Setup 

The holographic model that dualizes to a field theory in the presence of a helical lattice has 
the action [26] 


-^heiix = / 


R + 12- 


- - B AF AW. 


( 2 . 1 ) 


Here 5 (^ 1 / is the metric of a 5-dimensional asymptotically anti-de-Sitter spacetime including 
the 3 -|- 1 field theory dimensions and the additional radial coordinate r. R is the Ricci 
scalar of this metric. There are two field strengths: F^i, = d^Ay — dyA^ is the Maxwell 
field which accounts for the U{1) charge dynamics. The additional massive Proca field 
B^ generates the ‘helix U(l)’ with field strength W^y = d^By — dyB^, and supports the 
helical structure. In addition, there is a Chern-Simons term which couples the fields A^ 
and B^ with coupling constant n. In the above action, the AdS radius L has been set to 
one. Furthermore, Newton’s constant has been fixed to = 1/2. This can be achieved by 
redefining the remaining couplings such that becomes a total factor multiplying the 

action. 


To encode the 17(1) order parameter, we add to this action a scalar field with charge q and 
mass nip minimally coupled to 


•S'total = -S'helixT / d^+^x 


- |(9p-igApp 


The equations of motion following from the action (2.2) are 




( 2 . 2 ) 


(2.3) 


where 


- \9,uF\ 

Tfy^ = \Wp^Wy^ - ^gpyW^ - "^BpB, 


TjfJ = Re [{Vpp* +iqApp*){Vyp-iqAyp)\ - ]^gpy{\dp - iqAp\^ + ml\p\^), (2.4) 

are the energy-momentum tensors of the two vector fields A and B, and of the complex 
scalar p. Furthermore, we have the scalar equation 

0 = [(V^ - iqA>^){Vp - iqAp) - m^] p, 
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(2.5) 







Figure 2. Plot of the one-form uj 2 along the x-axis for one period. Being periodic with period 
2tt/p, lo 2 is not translationally invariant for p ^ 0. The vector field B = w{r)uj 2 acts as a source for 
the helix and imprints the helical, translational symmetry breaking structure on the system. 

and the Maxwell equations 

= iq [p*{d-' - iqA^)p - p{d^ + iqA’^)p*] + (2.6) 

[2d^{B^F^p) - . (2.7) 

Here is the totally antisymmetric Levi-Civita symbol in 5 dimensions with g0i234 _ 

As in [26], the wedge product in the action (2.2) is normalized such that the Chern-Simons 
term evaluated on the chosen Ansatz equals Fes = f drpKtv^a'/2. 

We now construct solutions to the equations that have the following properties. First we 
aim to study the system with the helix structure in order to break translational symmetry. 
For this purpose, the one-forms 

wi = dx, 

UJ 2 = cos{px) dy — sm(px) dz, 

LOs = sin(px) dy -t- cos{px) dz, ( 2 . 8 ) 

are introduced. They provide a basis for the spatial (x, y, z) part of the metric and the two 
vector helds and In Figure 2, one period of u }2 is plotted along the x-coordinate. 
The forms uj 2 and U 3 have the structure of a helix with periodicity In the following, we 
focus on the case m = 0, i.e. we are considering a massless helix held B. In our setup, the 
role of B is to introduce a lattice in a phenomenological way and thus break translational 
symmetry. Since this can be achieved with a massless helix held, m = 0 is chosen for 
simplicity. This choice follows [26]. Using these one-forms we make the Ansatz for the helix 
held B = Bfj_dx^ to be 
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B = w{r)uj2, 


w{oo) = A, 


(2.9) 




where r = oo denotes the boundary of the asymptotically anti-de-Sitter space. Since this 
Ansatz shows that By and do not vanish at the boundary, the held theory interpretation 
is that we explicitly introduce a source A for the operator dual to B, i.e. we are deforming 
the homogeneous theory by a lattice operator. A can be interpreted as the lattice strength. 
The held B extends along U 2 and therefore breaks translational symmetry in the x-direction 
for p 7 ^ 0. Via backreaction on the metric, this helical structure is imprinted on the whole 
gravitational system. This is manifested in a metric Ansatz [11] 

ds^ = -U{r) dt^ ul + e^^^M ^2 ^ ^2v3{r) ^2.10) 

From a technical point of view, the usefulness of this Ansatz is that it is compatible with 
the so-called Bianchi VIIq symmetry and is therefore guaranteed to be self-consistent — 
naively the spatial dependence in the one-forms should imply that all the components of 
the metric become spatially dependent. Thanks to the symmetry this is not so. Instead, 
all x-dependence is carried by the one-forms of Eq. (2.8) such that the resulting equations 
of motion are ordinary differential equations in the radial coordinate r. It is therefore con¬ 
sistent to assume that all helds are functions of r only, even though translational symmetry 
is broken. 

In the dual held theory the blackening factor U encodes the energy density and the Vi are 
related to the pressures in the system. In the hnite temperature phase, the metric function 
U has a zero at a hnite value of r, which dehnes the thermal horizon radius 

U{rh) = 0. ( 2 . 11 ) 

We will consider solutions which, for large values of r, satisfy 

U{r) = r'^ + ..., Xi(r) = ln(r)-|-..., fori = 1,2,3. (2.12) 

This guarantees that at the boundary (for r —>■ oo) the metric is of anti-de Sitter form, 
ds^ = dr^/r^ -|-r^(—dt^-|- dx^). Field theoretically, this means that the theory has an 
ultraviolet hxed point with conformal symmetry. The introduction of the helix through the 
source for the operator dual to By deforms away from this UV hxed point, as we will see 
below, c.f. Eq. (2.37). 

To model a held theory at hnite density, we make the additional Ansatz 

A = a{r)dt, a(oo) = p. (2.13) 

Field theoretically the operator dual to Ay is a conserved current j^. As this Ansatz again 
does not vanish at the AdS boundary it means that we introduce a source for the zero com¬ 
ponent of j^. We work in the grand canonical ensemble with chemical potential /i rather 
then a system at hxed density. 

We furthermore choose the charged scalar, which is the gravitational encoding of our su¬ 
perconducting order parameter, to respect the Bianchi VIIq symmetry, i.e. to be invariant 
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under the vector fields dual to the Blanch! VIIq one-forms in Eq. (2.8), 

uj'^d^p{t,x,y,z,r) = 0, i = 1,2,3. (2-14) 

Since the Bianchi one-forms (2.8) are linearly independent, this restricts the charged scalar 
to be at most of form p{t,r). For the background we choose a time-independent Ansatz 
p{r). Note that choosing a superconducting order parameter compatible with the helix 
symmetries corresponds, in condensed matter terminology, to the statement that the su¬ 
perconducting order parameter must respect the crystal symmetry of the underlying lattice. 
The resulting coupled ordinary differential equations for the functions a, w, U, vi, V 2 , V 3 
and p can be found in Appendix A. Making use of the U{1) symmetry associated with A^, 
p is chosen to be real. The total differential order of the equations of motion is 13, corre¬ 
sponding to the three second order equations for a, w and p, three second order Einstein 
equations and one first order Einstein equation, the constraint equation. 


It will be shown that this system exhibits a second order phase transition at finite tem¬ 
perature: above some critical temperature Tc, the system is in a state with a vanishing 
scalar field; below Tc a phase with a non-vanishing scalar field is thermodynamically pre¬ 
ferred. The order parameter of the phase transition is the vacuum expectation value {Op) 
of the operator dual to p, and the low temperature phase is characterized by breaking the 
global U{ 1 ) under which Op is charged. Since we will only consider the case of Op not 
being sourced explicitly, the U{1) symmetry is broken spontaneously. The superconducting 
phase transition can be understood in terms of an effective scalar mass becoming sufficiently 
negative: from the scalar field equation 


0 = p" + p' 


IT 

U 


+ v'l + v '2 + v'^ \ + p 


aV 

C/2 



(2.15) 


we can deduce the effective r-dependent scalar mass mlg = rUp — (f'a?/U. The second 
term shifts the effective mass squared towards smaller values and causes an instability as 
rUgjj becomes sufficiently negative. Both larger values of q and more negative values of 
favor the appearance of the instability. Most studies on holographic superconductors with 
a scalar condensate focus on tachyonic scalar masses corresponding to relevant operators.^ 
From studies on translationally invariant s-wave superconductors it is however known that 
condensation to a superconducting state also appears for a massless (and even irrelevant) 
scalar [39]. As in the case of translationally invariant s-wave superconductors studied in 
[40], it is expected that the zero temperature infrared geometry depends on the mass of 
the scalar field and that the massless scalar represents one of several distinct cases. In the 
following, we focus on the case nip = 0 and leave the case of general scalar masses to further 
studies. One reason for this choice will become clear in Section 4: We have only been able 
to construct extremal (zero temperature) geometries for this value of the scalar mass. The 

^Tachyonic scalar masses are allowed in asymptotically anti-de Sitter spaces as long as they are above 
the Breitenlohner-Freedman bound [38] > —d^/4. For such fields the tendency towards an instability 

driven by the negative mass squared is stabilized by the curvature of the AdS space. 
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mass of the scalar field nip is related to the scaling dimension A of the operator dual to 
p via A = 2 + y^ 4 + Therefore, the case nip = 0 corresponds to a marginal operator 
with A = 4. 

The Chern-Simons term couples the Maxwell field Ap and the helix field Bp. The authors 
of [26] point out that in their case the Chern-Simons coupling allows for the existence of 
a cohesive phase, i.e. a phase without electric flux through the horizon. This requires 
some field or coupling that sources the electric field outside of the horizon. The insulating 
geometry of [26] is of this type and it is stabilized by the Chern-Simons interaction, inducing 
the necessary charge density via the coupling to the helix field. In particular, for larger 
values of the Chern-Simons coupling k > 0.57 [26] the phase diagram simplifies due to 
the absence of an under RG flow unstable bifurcation critical point. Due to this added 
simplicity, and in order to be compatible with the results of [26], we focus in this work 
on the case of relatively large Chern-Simons coupling and choose k = YC2. Furthermore, 
in order to be able to compare the physics of our model better to the usual holographic 
superconductor [41, 42], we also investigate the case of vanishing Chern-Simons coupling 

K = 0. 

2.1.1 Asymptotic Expansions 

As an important step towards solving the equations of motion, asymptotic expansions are 
calculated at the thermal horizon rh and at the boundary, i.e. for large values of r. These 
are the two points where the boundary conditions on the fields are imposed. In addition to 
its defining property as the zero of the metric function U, the horizon is also a zero of 
the Maxwell potential a due to regularity {c.f. for example [43]) 

U{rh) = 0, a{rh) = 0. (2.16) 

The remaining fields are finite at the horizon. At the boundary r —)> oo, the conditions 

a = p, w = X, P = 0, tJj = ln(r), (2-17) 

are imposed. The first three conditions determine the sources of the operators dual to 
a, w and p. For general scalar field masses nip, the leading power of p at the boundary 
is where A is the scaling dimension of p. The massless scalar considered here 

has A = 4, such that it is constant to leading order in i/r. Therefore, the source of the 
operator dual to p is chosen to vanish by imposing p = 0 at the boundary. A solution with 
a non-vanishing scalar field in the bulk then breaks the U{i) symmetry associated with 
the Maxwell field spontaneously. As discussed before, the conditions on Vi ensure that the 
metric is asymptotically anti-de Sitter. For an asymptotically AdS space it is also required 
that U{r) = at the boundary. This condition, however, follows from the equations 
of motion (in particular the constraint, the sixth line in (A.l)) and does not need to be 
imposed explicitly. In order to determine the asymptotic horizon expansion respecting the 
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conditions stated above, we make the Ansatz 

a = Oi (r - rh) + a^ir - ruf ^-, w = Wq + (r - rn) ^ -, 

U = U^{r-rh) + U!^{r-rhf+ ■■■ , + '' • > (2-18) 

P = PO + Pli^" ~ ''"h) + ■ ■ ■ ) 

and solve the equations of motion order by order in r — rh- The horizon expansion has 
seven free parameters which are chosen to be 

iai,Wo^Po,Ui,v^^Q^). (2.19) 

All higher order expansion coefficients can be expressed in terms of these parameters. 
is related to the Hawking temperature by 

^ U'jrh) 

dvr dTT 

The asymptotic horizon expansion will be used to define initial conditions for the numerical 
integration of the equations of motion. 


At the boundary, a double expansion in i/r and in is carried out. The ln(r)-terms 

introduce a scale and indicate the presence of a scaling anomaly, which is related to the 
lattice structure. This will be discussed in more detail in Section 2.1.2, where the stress- 
energy tensor of the system is calculated. The following Ansatz is used as a building block 
for the asymptotic expansion at the boundary 


N N-k 

fNir) “fcb 

k=0 j=0 



1 Inr 

OOO + ( “10“ + ^ 01 — —I- 


1 


Inr 


+ 020^ + aii^r + ®02 


(Inr) 






( 2 . 21 ) 


The constant N defines the order of the expansion. Since both i/r and are small for 

large values of r, higher powers in the expansion are truly subleading. /Ar(r) is used to 
define an Ansatz for the matter fields a,w and p. The Ansatz for U is given by r^/ 7 v(r), 
and the Ansatz for the metric functions Vi by ln(r) -|- /Ar(r). This accounts for the fact that 
r^ is the leading behavior of U and of e^'“* for large values of r. Solving the equations of 
motion order by order in i/r, we obtain an asymptotic expansion with the leading terms 


d — p — 2 T • • • , 

2 e/3 + p^X‘^ln{r)/6 

^ ^ j^2 ' * * * ? 

W N , 52-p^A2ln(r)/12 
V 2 = ln(r) r-—— H-, 

^ n _ q^p^Pb 
^ r^ 12r® 


, /3 — p^Aln(r)/2 

w = \ + -— ^2 + ••• ’ 

w ^ , 51+P^A2ln(r)/24 

vi = ln(r) H- - -^- 

W ^ , 53+P^A2ln(r)/24 

113 = ln(r) -k -- 4 -— 4- 


( 2 . 22 ) 
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Here we find that 53 satisfies 


53 = -51 - 92- 


(2.23) 


As we will see later, the parameters gi are related to the pressure of the system and e is the 
energy density. The subleading mode of a denoted by u, is related to the charge density 
and /3, the subleading mode of w, to the vacuum expectation value of the operator dual to 
w. The leading mode of p has been set to zero and the subleading mode, p;,, is proportional 
to the vacuum expectation value {Op). In total there are 8 physical parameters at the 
boundary, namely 


A,/3,pfe). (2.24) 

There is another, non-physical parameter present in the asymptotic boundary expansion, 
which has been set to zero above to keep the expressions clear. This parameter is related 
to a shift in the r coordinate. It can be reinstated into the boundary expansion by the 
transformation 

+ 52-t52 + ^- (2.25) 

a gives rise to odd powers of i/r in the boundary expansion. These odd powers are in general 
present in the (numerical) solutions to the equations of motion. However, in contrast to 
the remaining parameters of the asymptotic expansion, a has no physical meaning. It 
can be thought of as an artifact of the coordinate choice. Note that the Ansatz of Eq. 
( 2 . 10 ), and therefore also the equations of motion, do not explicitly depend on the radial 
coordinate r. As a consequence, the system exhibits a shift-symmetry r ^ r + const, in 
the radial direction. The presence of a in the boundary expansion reflects precisely this 
shift-symmetry. 


2.1.2 Thermodynamics and the Conformal Anomaly 

The temperature of our strongly coupled field theory is given by the Hawking temperature 
of the bulk black hole, Eq. (2.20). The Bekenstein-Hawking entropy is calculated according 
to the area law as 


S 





Q^l{rh)+V2{rh)+V3(rh) y 


(2.26) 


Here denotes the area of the black hole horizon, jab is the induced metric and V = 
f d^x denotes the 3-dimensional volume of the field theory. In order to determine the 
grand canonical potential, we calculate the Euclidean (imaginary time) on-shell action. 
Applying appropriate integrations by part with respect to r, the Euclidean action reduces 
to a boundary term upon use of the equations of motion. We obtain the two versions 


.fbulk 

.fbulk 


V 
T 

V 
T 


-ae’'i+"2+^3 a' - -aKpw^ + y' 


r=ri, 


r=rh 


_yyi-V2+V3 ^^^V3+V2+V3 _^yyi+V2+V3 


-1 r=rj, 


J r=rh 


(2.27) 
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which will give rise to two expressions for the grand canonical potential n. A regularizing 
ultraviolet cutoff has been introduced. Eventually, after determining appropriate 

counterterms, the limit —?■ oo will be taken. The renormalized Euclidean on-shell action 
is given by Iqs = /bulk + /gh + let-, where /qh is the Gibbons-Hawking boundary term and 
let stands for the counterterms necessary to make Iqs finite.^ The factor i/t = (5 originates 
from the r-integration and V = f d^x is the three-dimensional volume. Note that the 
second expression for /bulk receives no contribution from the horizon since U{rh) = 0. The 
Gibbons-Hawking term is evaluated as 


/gh = -2 


drd^x^ 


V 

T 


gUl-|-U2-|-U3 jjf _ g'^’1+^’2+^’3 


v[ - 2 //e^i+’' 2 +t ;3 _ 2 [/e^i+^ 2 +i <3 


(2.28) 


Here 'jab is the induced metric at r = r^, and n = \/Udr is the outward pointing normal 
vector to the surface at r = r^. The sum /bulk + /gh is divergent in the limit r;, —)• oo. 
Using the asymptotic expansion of Section 2.1.1, the diverging terms can be written as^ 


/bulk + /gh 



In (rb) 


1 

2 


pKfiX^ “t" 2e -{- 2/rz2 -j- O 



(2.29) 


Here the first expression for /bulk of Eq. (2.27) was used. The divergent terms are, however, 
the same for both expressions. They only differ in terms that are finite at the boundary. 
There are two types of divergences: the — 6 r^-term is due to the infinite volume of the 
asymptotically anti-de Sitter space and can be canceled by a counterterm proportional 
to the volume f drd^x ^/7 of the surface at r = r/,. The second type of divergence is 
logarithmic in and, in order to cancel it, a counterterm proportional to In(rb) needs 
to be introduced. As shown later, this causes a scaling anomaly. Roughly speaking, the 
logarithm present in the counterterm introduces a scale and therefore breaks conformal 
symmetry. The total on-shell action Iqs = /bulk + /gh + /ct can be made finite by means of 
the counterterm 


/ct = 


= J drd^x-y /7 


6 - - ln(rb)lT“^iy,b 


(2.30) 


In this expression, indices are contracted using the induced metric ■jab- Since the scalar 
field p does not cause any new divergences, these counterterms are identical to the one of 
[26].^ 


According to the AdS/GFT correspondence, the field theory partition function is related 
to the on-shell action by Z = exp(—/ qs). Therefore, the grand canonical potential can be 
expressed as H = T/qs. Using the asymptotic boundary expansion of Eq. (2.22), H can be 

^For a review of holographic renormalization, see [44] . 

®For the sake of clarity, the shift parameter a is set to zero in this expression. It can be reinstated using 
the transformation of Eq. (2.25). 

“^In the corresponding expression of [26], an additional term proportional to ^7 ln(ri,)F“^F’af, is included. 
This term is present for general gauge fields but vanishes if has only a time component. 
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written as 


n 

V 

n 

V 


€ + 2(11^ — -kX'^^p — 

4 

CX 6 

451 + ----/3A- 
1 d 6 


Q'm(rh)+V2irh)+V3(rh) JJ> ^ 


(2.31) 


corresponding to the two expressions for /bulk of Eq. (2.27). The non-physical shift- 
parameter a, which is explicitly indicated here, can be removed by the redehnition (2.25). 
The hrst expression for 1/ can be brought into the more familiar form 


5/ 

V 


= e — pn — Ts, 


(2.32) 


by rewriting the horizon contribution in terms of the temperature T and the entropy S and 
by dehning the particle density as 


n = -2v-\ -(2.33) 

where s = ^/v denotes the entropy density, and e is the energy density. This will be 
conhrmed by an explicit computation of the held theory stress-energy tensor. Since the 
scalar p has a vanishing source, there is no explicit scalar contribution to the grand canonical 
free energy.^ Therefore the above expression agrees with the one given in [26]. 


The expectation value of the held theory stress-energy tensor [45, 46] is calculated from the 
extrinsic curvature Bab at r = rb, its trace 9 = ^-nd the real-time counterterm action 

Set U-S 

Using the asymptotic boundary expansion, the resulting expression for {Tab) can be ex¬ 
pressed in terms of boundary parameters as 


with 


{T) = 


(e 0 0 0 \ 

0 8^1 - 4AV + i 0 0 

f) n T T 

u u ^yy 

\0 0 Ty, 


'^yy = \- ^9i - ^AV + {^9i + 892 + cos( 2 px), 

Tyz = - ^ 4 c/i -/ 852 + sin(2px), 

^ ^ \ 2 J 2 I 1 \2 2 


Tzz = -- 451 - -X p - ( 4c/i + 8 g 2 + -X p ] cos(2px). 


(2.35) 


(2.36) 


® The grand canonical free energy does, however, implicitly depend on the scalar. When existent, a 
solution with a non-vanishing scalar has lower free energy than the normal phase solution, as the superfluid 
transition is second order in our model. 
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For p ^ 0, the stress-energy tensor is spatially modulated and anisotropic. Indeed, e = Too 
is the energy density, and the parameters gi are related to the pressure terms Tn, T 22 , 
and r 33 . The trace of the stress-energy tensor {T°‘a) — i'^ab) determines the conformal 
anomaly. It is given by 

(n> = - 2 ^V. (2.37) 

In the presence of the lattice, i.e. for p ^ 0 and A 7 ^ 0, the conformal symmetry of the 
ultraviolet fixed point is broken by an anomaly. The fact that (T“^) is proportional to has 
the intuitive interpretation that the lattice momentum p introduces a scale and therefore 
breaks conformal symmetry.® In addition to the temperature T and the sources p and A, 
the system is characterized by the dimensionful scale p. Out of these four dimensionful 
parameters, we can form the dimensionless ratios Ym and p/im. These ratios will be used 
to parameterize the system (in addition to the scalar charge q, which is dimensionless). 
Note however that p itself is a parameter in the state of the boundary field theory, but 
not a source switched on and hence, as explained in detail in Appendix D, should not be 
counted as an independent UV integration constant. 


2.2 Phase Transition with a Scalar Order Parameter 


Numerically, we find that above some critical temperature Tc, there is only a solution to the 
equations of motion with a vanishing scalar field. ^ For T < Tc, a second branch of black hole 
solutions arises which has a non-vanishing scalar field and therefore a non-vanishing vacuum 
expectation value {Op). Since the source of Op is set to zero, the condensate {Op) breaks 
the U{1) symmetry associated with the Maxwell field Ap spontaneously. By comparing 
the free energy in the grand canonical ensemble, we find that the solution with a scalar 
condensate is thermodynamically preferred, c.f. Figure 3. The phase transition is second 
order, since the derivative of the entropy is discontinuous at the transition temperature. 
The order parameter of the phase transition is the vacuum expectation value {Op), which 
is proportional to the boundary mode pb of p. It is plotted as a function of the temperature 
in Figure 3. The numerical data are consistent with a mean field behavior of the order 
parameter near the transition temperature, {Op) ~ This is expected in the 

large-limit, which is intrinsic to our holographic model. For each set of parameters 
{p/ 11 ,^/ii,q), the transition temperature can be determined as the temperature where the 
order parameter {Op) vanishes. In this way, the phase diagram of the strongly coupled field 
theory can be studied numerically. The transition temperature increases as a function of q 
as is shown in Figure 4. This can be understood in terms of the effective scalar mass 


"leff = - 


U 


(2.38) 


As the scalar charge increases, the effective mass squared becomes more negative which fa¬ 
vors the instability. Therefore, the phase transition already appears at higher temperatures. 

®Note that the transforms under the scaling transformation (D.l) with weight 4, as it should. 

^The numerical method used is a simple shooting method adjusting the free parameters of the horizon 
expansion in such a way that we arrive for certain fixed values on the boundary. For a detailed exposition, 
see Appendix D. 
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Figure 3. In the left panel, the free energy for solutions with and without a scalar condensate are 
shown. The blue curve corresponds to solutions without and the red curve with a scalar condensate, 
respectively. Having a lower value of the grand canonical potential fl, the solution with a scalar 
condensate is thermodynamically preferred. The phase transition is second order since the derivative 
of the entropy is discontinuous at T = T^, see inset. In both graphs k = 9 = 5, vjp, = 5.2, and 

Ym = 3.3. In the right panel, the scalar condensate order parameter (O) is shown for p/m = 4.8, 
5.0, 5.2 and Ym = 3.0, 3.3, 3.0, 3.9. The graph shows that the strength of the order parameter is 
affected much stronger by changes in Ym than by changes in p//^: different A leads to each of the 
four curves shown, while the different vj^i lead to nearly the same curves for each value of Ym- 


K = 0 @ p/m = 1.2 K = 1/V2 @ p/m = 1.2 




Q q 

Figure 4. The superfluid phase diagram as a function of the scalar charge q close to the minimal 
Tc-value in the phase diagram, i.e. p/^i = 1.2 for k = 0 (left figure) and k = (right figure). 
As one can see the critical temperature is decreasing with decreasing scalar charge q or increasing 
backreaction ~ i/q, respectively. This result agrees nicely with the plain s-wave superconductor 
behavior. Moreover for k = i/Vz and for sufficiently high values of Ym (color coding of the lines 
represent the following values of Ym = !> 2, 3, 4, 5, 6 , 7, 8, 9, 10) we expect a quantum critical point, 
where the superconducting phase transition happens at zero temperature. On the other hand, for 
K = 0 it seems that the curves are asymptotically approaching zero charge as Y = 0, and thus 
there may be no quantum critical point in the q direction. 


The same behavior was found for translationally invariant holographic s-wave superconduc¬ 
tors in [42]. 

Analyzing the transition temperature as a function of P/^i and Ym for different q reveals a 
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K = 0 q = 6 



Figure 5. The superfluid phase diagram computed for k = 0 and scalar charge (or inverse back- 
reaction) <7 = 6 in the top panel and for k = Y %/2 and g = 5,10 in the lower panel. Below the 
respective Tc-curve, the condensate is non-zero indicating a superfluid phase. Different colored 
curves correspond to different values of ^/^l with the following color coding: 0, 0.3, 0.6, 0.9, 1.2, 
1.5, 1.8, 2.1, 2.4, 2.7, 3, 3.3, 3.6, 3.9, 4.2, 4.5, 4.8, 5.1, 5.4, 5.7, 6. Note that the case of k = 0 the 
minimum of the critical temperature is moving from p/fj. = 1.8 at Ym = 1-2 over p/a* = 2 at Ym = 3 
to p/fi = 2.4 at Ya* = 6 > whereas in the case of k = Y %/2 the minimum is fixed respect to A. However 
for different values of the scalar charge g = 4,..., 10 we And a slight increase starting from p/h = 1.2 
and ending at p//j = 1 -6- 


more interesting structure as shown for q = 5 and q = 10 in Figure 5. In the case of k = 0 
the critical temperature is observed to decrease monotonously as a function of Ym- As a 
function of p/m; it first decreases for small pf/i, then assumes a minimum at p//i ~ 1.7, which 
is slightly shifted towards larger values i.e. p/ti ~ 2.3 for increasing Yp Y 10, and then 
returns again to the homogeneous value Tc{p/p. = 0) for large values of p/^l. This minimum 
is more pronounced for larger values of Ym which shows that larger values for the source of 
the helix held indeed increases the effect of the lattice on the system, while the helix mo¬ 
mentum dependence has a smaller effect on the transition temperature. This is consistent 
with the expectation that generally, it is the depth of a lattice of potential valleys which 
influences the physical behaviour more than the lattice constant or spacing between the 
individual potential depths. However, for large values of p/^i the critical temperature seems 
to, at least for k = 0 , asymptotically approach the p = 0 value which might imply that 
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Tc{p oo) ~ Tc{p = 0). 


In the case of k = there is not only a minimal value of Tc for p/ti ^ 1.2 that is robust 
under changes in Y/i, but allows for higher values of Tc above p/ii ~ 2.6,..., 3, as displayed 
in the lower row of Figure 5. Curiously, there the behavior of ^//i is inverted, i.e. with 
increasing Y/x the transition temperature, Tc is increasing. However, this regime might very 
well lie far outside of the range of applicability of condensed matter physics since p > p, 
suggesting that the “energy stored” in the lattice exceeds the chemical potential. It appears 
that the critical temperature is unbounded from above unlike in the k = 0 case where it 
seems to be bounded by its initial value at p = 0 . 

Finally, note that in both cases ac = 0 and k = 1 /V 2 the data suggests the existence of a 
quantum critical point for high values of Ym > 6 and strong backreaction, i.e. q < 5, similar 
to the observations made in [30] . One may hence speculate that for a finite range of q there 
exists a critical value of Y/^ zero temperature where the superfluid phase breaks down 
above a certain p/ii. We will analyse this possibility further in future work [37]. 

3 Optical Conductivity 

A main focus of our work is to study the optical conductivity in the x-direction, in which 
translation symmetry is broken for p ^ 0. The conductivity is given by the Kubo formula 
in terms of the retarded Green function of the current operator, 

axiio) = hm-. (3.1) 

k-!>0 lOJ 

3.1 Numerical computation of the optical conductivity 

We calculate the retarded Green function using the well-established extension of the gauge/- 
gravity correspondence to real-time problems pioneered in [47]. Accordingly, the conductiv¬ 
ity is determined by linearized perturbations around the background solutions. In a physical 
picture, the system is slightly perturbed around thermal equilibrium by an external force, 
the external electric field. The conductivity describes the induced response of the system 
to these small perturbations. Since the gauge field corresponds to a conserved current 
on the field theory side, the electric conductivity is related to the perturbations of A^. In 
particular, the fluctuation 5A = A dx governs the conductivity in the x-direction. There 
are certain fields coupling to 6 A = Adx. To determine them, we write all fields including 
the metric as a sum of the background solution and a perturbation. Then, the action of 
(2.2) is expanded to second order in the perturbations. In this way, a quadratic action Sq 
for the perturbations is obtained, which determines the linearized equations of motion for 
the perturbations. By analyzing Sq, all fields coupling to SA = Adx are determined.® The 

®The full set of the most general couplings at vanishing momentum is shown in the Appendix A.l, 
Table 1. 
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block of coupled perturbations containing 5A = Adx is 


5 (ds^) = + h 2 'i{t,r)uj 2 0 ws + hri{t,r) dr ® 0 Ji, 

5A = A{t,r)ui, (3-2) 

6 B = B{t, r)iOs. 


The remaining perturbations can be set to zero, consistently. The fluctuation fields are 
chosen to depend on r and t only since the conductivity is evaluated in the limit of vanishing 
spatial momentum for the Kubo formula (3.1). The equations of motion for the above 
perturbations are obtained by varying the quadratic action. After variation, we impose a 
radial gauge in which hri = 0. In this gauge, the equation for hri becomes a constraint for 
the remaining fields. Furthermore, the equations become ordinary differential equations in 
the radial coordinate after Fourier transforming the time coordinate. In total, we obtain 
one first order equation (the constraint originating from the hri-equation after choosing 
radial gauge), and four second order equations for A,B,hti and / 123 . One of the second 
order equations can be replaced by the constraint, hence the total differential order of the 
system is 7 = 1 -|- 2 x 3. The equation of motion are given in Appendix A.l. The asymptotic 
expansions of the fluctuation fields, which are necessary for obtaining numerical solutions, 
are discussed in Appendix B. As usual we implement inf ailing wave boundary conditions 
at the thermal horizon r^- For A, we find that 

A = (r- (aS + Al{r - r^) + ■ ■ ■) . (3.3) 

The exponent in the prefactor can assume the values which correspond to out¬ 

going and infalling waves, respectively. This can be seen by taking into account the phase 
factor of the Fourier transform. 


(r - 


giaj(±r-t) 


with 


ln(r - rfe) 

(dvrr) 


(3.4) 


The infalling solution, i.e. the one with the minus sign, should be used in order to obtain 
the retarded Green function. The leading behavior of A near the boundary is found to be 


A 


4 + 


A 2 + Aga;^ln(r )/2 


+ 


(3.5) 


For calculating the conductivity, we have to identify the degrees of freedom coupling to 
A by imposing gauge invariance. Even after imposing radial gauge, in which all radial 
fluctuations vanish, there are still residual gauge transformations left. These consist of 
the diffeomorphisms and U{1) transformations that do not change the radial gauge = 
Ar ^ Br ^ 0. The residual gauge transformations are worked out in detail in Appendix 
B.l, following a calculation carried out in [48] in the framework of the holographic p-wave 
system. We find that the relevant physical fields are (i) A, which is already gauge invariant, 
(ii) /i 23 , which is gauge invariant at the boundary r —>■ 00 , and (iii) the linear combination 


Q = —iuiB + wpe hti. 


(3.6) 
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The field hti is not gauge invariant and does therefore not carry physical degrees of freedom. 
In order to calculate the Green function corresponding to A, we impose the condition that 
the remaining physical fields, h 23 and Q, have no source term, i.e. that their leading modes 
for r —)• oo vanish.® In this case the renormalized on-shell action for the fluctuations (c./. 
Appendix E), expressed in terms of the asymptotic modes of A, is^® 


•S'on-; 


on-shell 


dw 


d3xA^(-u;) 




(3.7) 


In this expression, the frequency dependence of the modes is indicated explicitly. The Green 
function does not follow directly from (3.7). Following the prescription of [47], one needs 
to analytically continue the kernel in (3.7). It follows that 


G^,(u;,0) = 2 




(3.8) 


and, using the Kubo formula (3.1), 


a{uj) 


2^2(0;) \ 0J 
ia;Ag(a;) 2 


(3.9) 


The numerical steps in calculating the conductivity for a given solution to the background 
equations of motion are described in Appendix D. 


3.2 Comparison to the Drude-model and the Two-Fluid-model 

For holographic metallic systems in homogeneous translation invariant backgrounds, one 
finds an ideal metallic behavior related to the conservation of momentum. Thus, strictly 
speaking the Drude model is not applicable. In the presence of a lattice, momentum is not 
a conserved quantity and therefore charge carriers can dissipate their momentum within a 
typical timescale r by interactions with the lattice. According to the Drude model (c./. for 
example [36] for a review and more details). 


a{uj) 


<700 

1 — lUJT ’ 


Re 17 ( 0 ;) 


O'DC 

1 -I- (wr)^’ 


(3.10) 


the dissipation time scale r is inversely proportional to the width of Re((T) near ui = 0, and 
the Drude peak can be seen as a direct consequence of the translation symmetry breaking 
lattice. The same reasoning carries over to holographic superconductors. However, in the 
absence of a momentum dissipating mechanism such as a lattice, the holographic system 
describes an ideal metal in the normal phase and a mixture of a holographic superconductor 
and a remaining ideal metal in the condensed phase. In the limit of restored translational 
symmetry, the Drude peak degenerates into a delta peak at w = 0 In our helical setup, 

^Alternatively, we can make use of a method devised for treating holographic operator mixing, as ex¬ 
plained in Appendix E. 

^°Only the boundary contribution is indicated. According to the prescription of [47], the horizon contri¬ 
bution is to be discarded. 
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q = 6 tT = 1.22 


q = GkT = 0.41 7 ; 




Figure 6. The left panel displays the optical conductivity in the translationally invariant system, 
i.e. p/fj. = 0, or = 0, respectively. The curves correspond to q = 6 and T = 1.22 71,. The '^/oj 
pole indicating an infinite DC conductivity is also present in the normal phase due to momentum 
conservation. Accordingly, there is no momentum relaxation possible and hence there is no Drude 
peak at small frequencies. In the superconducting phase, shown in the right panel by the blue lines, 
the delta peak is enhanced by the superconducting degrees of freedom and the system is a mixture 
of an ideal metal and a superconductor for finite temperatures T = 0.41 Tc, as shown in the inset. 
Comparing the w —1 0 value of the red line in the left panel u}lu\a /« 0.8 to the black in 
the inset in the right panel wlmcr//r^|^_g « 0.5, the “strength” of the ideal metal is reduced but it 
remains present, whereas the superfluid strength obtains a value of approximately 1.9. 

translational symmetry can be restored by setting p = 0 and/or A = 0.^^ In this case, the 
helix field w decouples from the system and we obtain the classical holographic model of 
an s-wave superconductor as introduced in [41, 42]. We can understand the translationally 
invariant case as a limit in which the relaxation time r tends to infinity: for lot ^ 1, the 
Drude conductivity reduces to a i/iv pole in Im((T), 

.ctdc 
a = 1 -, 

LOT 

which indicates an infinite DC conductivity as explained below (3.12). Charge carriers which 
are accelerated by the external field cannot dissipate their momentum and therefore the 
resulting response is infinite. Consequently, there are two physical mechanisms leading to an 
infinite DC conductivity in the translationally invariant case. First, there is a contribution 
for T < Tc which is caused by superconductivity. Additionally, there is a contribution due 
to momentum conservation. It is, therefore, necessary to break translational symmetry 
in order to determine the superconducting degrees of freedom separately. In Figure 6, we 
consider the translationally invariant case, and indeed observe the absence of a Drude peak, 
and the presence of a delta peak S{lo) both in the normal as well as superconducting phases. 

^^Setting p = 0 restores translational symmetry but the system is placed in an external magnetic field. 
In order to fully restore the plain holographic s-wave superconductor, the helix field needs to vanish, i.e. 
A = 0. The magnetic field however points in the direction of the helix director (the x-direction), and hence 
does not lead to a gap in the conductivities considered in this paper. 


(3.11) 
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normal metallic phase 


condensed phase 




Figure 7. The optical conductivity for q = 0, k = 0, v/ii = 2.4/^ Yt* = 4 and T = 0.57 Y 
in the normal phase (left panel) and in the thermodynamically preferred condensed phase (right 
panel). The real part of the optical conductivity exhibits a broad Drude peak for small frequencies 
in the normal phase related to the non-conserved linear momentum in the Bianchi type VII helical 
background. Comparing the normal phase real part of the optical conductivity (blue line) in the 
left panel with the condensed phase real part of the optical conductivity in the right panel, the 
developement of a soft gap becomes apparent. The gap scales algebraically with the temperature. 


Turning on the helical structure p ^ 0 and A Y 0; linear momentum is no longer conserved^^ 
and we hnd, at last for a weak helix Ym <C 1, a bona hde Drude-model behavior. In 
Figure 7, the optical conductivity is shown for p/m = 2.4, Ym = Ij and q = 6 for a certain 
choice of temperature in the broken phase and for the transition temperatnre Tc- For 
small frequencies, a Drude peak in the real part of the conductivity is observed both in the 
normal phase at Tc and in the snperconducting phase. Since we are not strictly at zero 
temperature, there is a remaining small Drude-like peak after subtraction of the '^/uj pole. 
Figure 8 shows the small-frequency regime of the optical conductivity a and a corresponding 
ht to the Drude model. Furthermore, in the condensed phase the imaginary part of the 
optical conductivity exhibits a pole for T < Tc indicating a delta peak in the real part 
of the optical conductivity related to an inhnite DC conductivity, which is a characteristic 
of superconductivity. This can be inferred from the Kramers-Kronig relation 


T / \ 2a; _ 

Im(T(a;) =- V 

vr 


. Recj(a;) 
da; 


0)2 — a;2' 


(3.12) 


According to this relation, a ^/uj pole in the imaginary part of the conductivity is related to 
a delta function at zero frequency in the real part by 


Imcr(a;) = — y Iiea{uj) = psd{uj). (3.13) 

vr oj 

This equation dehnes the superfluid density ps as the coefficient of the zero frequency delta 
fnnction in Re((T).^^ 

^^The canonical momentum related to the Bianchi VII group translations is still conserved, albeit it is 
not accessible on the boundary field theory. 

^^This value of p/m allows for the minimal Tc in the case k = 0 c.f. Figure 5. 

^^In the conventions of [33, 34, 36], pa is defined via Refer) = Ps5{lo)/8, i.e. it differs by a factor of 8 from 
the definition used here. 
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Figure 8. The small-frequency behavior of the optical conductivity in the metallic phase for 
q = 6, K = 0, p/ll = 2.4 and Ym = 1- The solid lines are fits of the Drude model (3.10) to the 
numerically determined optical conductivity. The normal phase solution is given by ctuq = 13.6873 
and r = 17.4352. In the superconducting phase the i/w Pol® has been subtracted from Im cr and 
there is a remaining “residual” Drude-like peak cr''®® as shown in the left panel. As explained in 
the main text and in Figure 13, this residual contribution to the two-fluid model (3.14) shows the 
coexistence of the superconducting phase and a “normal” holographic metal. 

As shown in Figure 9, a small Drude peak remains present in the superconducting phase. To 
describe the system, it is thus necessary to apply the two-fluid model [4] , which supplements 
(3.13) with the metallic Drude model defined in (3.10), 

Recr(a;) = + Ps5{uj) = (xn{T) ^ -h '^XsiT)6{uj)] , (3.14) 

\ 1 -h 2 / 

where Xn{T) describes the Drude-like contribution resembling a normal fluid and Xs{T) the 
superconducting contribution. In the normal state, we have Xn{T > Tc) = and XsiT > 
Tc) = 0, whereas a pure superconducting state would be described by Xn{T < Tfl) = 0 and 
Xs{T < Tc) = Ps-^^ Due to charge conservation Xn{T) + Xs(T) = n.^® 

Moreover, from Figure 9 we observe that the conductivity in the superconducting state 
develops a gap at low frequencies, i.e. Re(iT) is significantly suppressed. This gap is a 
characteristic of a superconducting system; it indicates that low-energy charged degrees of 
freedom have condensed into the delta function at cj = 0. An important issue is whether 
in (3.14) vanishes in the limit T —?■ 0 for frequencies below Wgap- This would imply 
Us = Ps at T = 0, with Hs the thermodynamic density. In general, the following scenarios 
are possible. One possibility is the presence of a hard gap, in which at low frequencies 
0 < w < cjgap we have an exponential suppression a{u},T) ~ exp((a; — ujga_p)/T). On the 
other hand, for a soft gap there is an algebraic (power law) scaling a{uj,T) ~ In 

^®In order to restore the proper units of the two-fluid model, the charge density is given in units of ®7-*. 
i.e. the number density and the charge density are related by Ucharge = e^nnun,ber/»ri*. Note that we work 
with charge densities and not number densities throughout the paper, as the quantities e and m* are not 
directly accessible in holographic models. Furthermore, this choice of dimensions has the advantage that 
the superfluid strength ps and the charge density n have the same units (in natural units). 

^®Throughout the paper n denotes a general charge density, while n„ denotes the charge density in the 
normal phase, and Ua the charge density in the superfluid phase. 


23 - 











Figure 9. The left plot shows the “insulating” low-temperature phase optical conductivity with 
vanishing imaginary part as w —>■ 0, which indicates the absence of a delta peak at w = 0 and 
allows for a proper Drude-like shape. In the right plot the superfluid phase shifts “spectral weight” 
from the finite frequency part which condensed in a a; = 0 delta peak, thus generating the gap in 
the superconducting state in the region “//x G [0,1.5]. Note that this is not a hard gap due to its 
algebraic temperature scaling on approaching the zero temperature limit and that there is a non-zero 
contribution coming from a residual Drude-like peak tr''®® even at small T in the superfluid phase. As 
expected from the Kramers-Kronig relations the imaginary part of the optical conductivity displays 
a i/tj pole. Both plots are generated for k = ^/V 2 , q = 6, p/m = 1.2, Ya* = 1-9, and T = 0.46 Tc. 


this case, it is much harder to determine numerically whether there exists an additional 
constant contribution cto to a{uj,T). As we discuss below, for the model considered in this 
paper, we find an algebraic scaling. Moreover, by calculating and ps independently, we 
find that Ug = Ps for T —)• 0 to good numerical accuracy, at least for small p/^l and Ym- This 
implies that (Tq = 0 in this case, i.e. our system exhibits a soft gap. 

For translationally invariant holographic s-wave superconductors in 2 -|- 1 dimensions, it is 
known that, even though highly suppressed, Re(cj) remains nonzero at T = 0 but finite fre¬ 
quency [40] and one finds an algebraic scaling. On the other hand, p-wave superconductors 
are reported to exhibit a hard gap [49]. Also for the helical lattice, by straightforwardly 
generalising the low frequency analysis in the Appendix of [26], we conclude that the gap 
scales algebraically in oj. 

In addition, we compute ps and Ug individually at very low temperatures and find indica¬ 
tions that they agree to good numerical accuracy, ps is read off from the zero frequency 
pole of the conductivity, while the thermodynamic density Ug is obtained from (2.33). In 
Figure 10, the superfluid density pg and the charge density Ug, c.f. Eq. (2.32), are plotted 
as a function of temperature for two sets of parameters. The superfluid density, being a 
measure for the superconducting degrees of freedom, increases as the temperature is lowered 
beyond Tc and it vanishes for T > Tc- Of course, in order to finally conclude whether ps 

^^Normalization of the charge density entails a rescaling by compared to the usual number den¬ 

sities in the Drude model, thus the dimension of ris is (length)”^ matching the dimension of ps = 
Y 2 lim,j_>o Im namely \ps\ = (length)”^. 
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K = 0, g = 6, p/fi = 3 & Y/i = 0.3 


K = 0, g = 6, p/m = 3 & = 1.5 




K = 9 = 5, p/p = 0.2 & Vm = 0.3 


K = i/v/ 2 , g = 5, p/p = 0.2 & Vp = 1-2 




Figure 10. The superfluid density ps and the normalized^^ charge density n as a function of 
For T > Tc, the superfluid density vanishes, i.e. the DC conductivity is finite. As T is lowered, 
the superfluid density p^ increases, similar to the order parameter, c.f. Figure 3, and curiously the 
normalized charge density rig as well. The gap between p^ and rig is independent of the value of 
p/ll., but increases with increasing Ym- For small ^/ii, it is suggestive that at T — 0 the superfluid 
density and the condensed phase charge density coincide. Thus, the longitudinal response i.e. the 
plasma frequency in the superconducting phase is sufficient to determine the superfluid strength. 


agrees with Ug for T = 0, as the extrapolation of our data suggests, we need to carefully an¬ 
alyze the zero temperature transport properties, which we are planning to do in future work 
[37]. Nonetheless, for small p/m Vm the difference becomes sufficiently small already at 
finite temperatures about T ~ 0.6 Tc. The difference between ps and Ug at this temperature 
seems to be independent of the helix pitch^®, parametrizing the helical lattice constant, but 
grows with increasing Y/x. This difference may be accounted for by the residual contribution 
in the condensed phase for a; —)• 0, which is not added to the zero mode delta peak. The 
optical conductivity of high superconductors is known [50] to feature residual absorption 
at very small frequencies, which gives rise to an additional contribution to the imaginary 
part of a{u,T < Tc). For small helix strengths Ym 1, this residual part can be read off 
by a simple Drude-fit inside the superconducting gap as shown in Figure 8. The spectral 

^®Technical problems arise for p/m ^ 1 due to numerical instabilities. 
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Figure 11. For k = 0, v|^^ = 0.8 and Ym = 3 we find that the “superconducting gap” is replaced by 
a plateau-like decay in stark contrast to the right panel of Figure 9, which shows the typical form 
of the optical conductivity. It is suggestive to compare this result to Figure 16 in [29] where the 
authors find a similar behavior for low temperatures. Again, the color coding for the temperature 
is as follows: = 0.47, 0.42, 0.37, 0.33, 0.28, 0.23, 0.19, 0.14, 0.09, 0.05. Counter-intuitively, the 

gap is drastically diminished with decreasing temperature, whereas one would expect the gap to 
become more pronounced for low temperatures. 


weight inside the residual Drude peak accounts exactly for the difference between ps and 
Us- On the other hand, for larger helix strengths, he. stronger momentum dissipation, the 
gap cannot be accounted for by the residual spectral weight inside the superconducting gap. 
We discuss two possible reasons for this behavior in Section 5.2. 

3.2.1 Intermediate and High Frequency Regimes 

Let us also discuss the behaviour of the optical conductivity in further frequency regimes. 
First we note that in the intermediary frequency regime T <C cu <C /i, we have not been 
able to find a scaling law of |cj(ti;)|. Such a scaling has been observed in the strange metallic 
phase of the cuprates and interpreted as a consequence of quantum criticality in [51]. While 
some holographic models [1, 4, 5, 52] seem to show such a scaling, others [26] do not, and 
our model seems to be in the latter class. So far a theoretical understanding of the origin 
of this scaling regime in holographic models is still missing. 

Concerning the nature of the superconducting gap, there are two more intriguing features: 
in the case of k = 0, we find in the vincinity of particular parameter values such as e.g. 
v/p, = 0.8, = 3, plateau-like solutions where the energy scale of the gap as a well as 

the superfluid density ps is drastically reduced, see Figure 11. Curiously, these solutions 
seem to arise for very low temperatures contrary to the intuition that the gap should grow 
with decreasing temperature, c.f. the orange line compared to the purple line in Figure 11. 
A similar behavior has been found in Figure 16 of [29], although in our case there is no 
Drude-like peak for higher temperatures close to the transition temperature, which may 
be attributed to the fact that almost all degrees of freedom have condensed.^® The same 

^^See also the discussion below Figure 13 and the FGT-Section 3.3. 
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K = 1/v^ & g = 6 


K = 1/^ & g = 6 




Figure 12. Different shapes of the optical conductivity with additional features compared to 
Figure 9 for p/fi = 0.2 and Ya* = 2. In the case of k = i/v^ for decreasing temperature the gap 
seems to get smaller and a new peak arises, interestingly at lower temperatures. The color coding 
for the temperature is as follows: Yt = 0.41, 0.37, 0.33, 0.29, 0.25, 0.2, 0.16, 0.12, 0.08, 0.04. Note 
that in all cases there is a delta peak in the real part indicating the superfluid phase. 

intriguing features are seen for k = y \/2 in the vicinity of particular parameter values 
such as the ones shown in Figure 12. Compared to the more generic result of the optical 
conductivity shown in Figure 8 and 9, the gap seems to develop a new feature in the 
intermediary frequency regime that resembles a peak at a non-zero frequency (as shown in 
Figure 12 at Ym ~ 1-78) at lower temperatures in the superconducting phase. Again not 
only the gap is reduced but also the superfluid density; this is evident since the Ferrell- 
Glover-Tinkham sum rule holds, as proven in the following Section 3.3. These curiosities 
appear to be due to the contributions of additional resonances below the gap. We plan to 
investigate these further in future work [37]. Finally, in all cases considered, we find that 
for large frequencies, w 3> /r, the real part of the conductivity is proportional to ui. This 
behavior is a property of the ultraviolet fixed point of the field theory: since a has energy 
dimension one, it scales linearly in u once the frequency is larger than any other scale of 
the system. 

3.3 Ferrell-Glover-Tinkham Sum Rule 

Sum rules are exact identities following from the analytic structure of Green functions. The 
Ferrell-Glover-Tinkham sum rule [53, 54] can be expressed as the integral over the real part 
of the optical conductivity, being a constant regardless of the details of the system, 

rojc 

/ dw Re cr(a;) = const., (3.15) 

Jo 

Note that the integral includes possible contributions from the lower bound in the form of 
a delta function at cj = 0. Furthermore, while actual condensed matter systems typically 
become transparent at frequencies larger than the typical electronic energy scales and hence 
the integral in (3.15) and (3.16) converges, in holographic setups it is generically UV di¬ 
vergent due to the UV conformal fixed point behavior of Re < 7 ( 0 ;). In order to regulate this 
divergence, we introduced a UV cutoff frequency ujc, to be taken larger than max(r,//). In 
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appropriate formulations of the sum rule such as (3.17) the UV divergence cancels between 
the two integrals, and the regulator can be removed.Physically, the sum rule expresses 
the conservation of charged degrees of freedom, which are measured by the spectral weight, 
i.e. the area under Re(iT). For example, in the normal phase, Eq. (3.15) allows to identify 
the plasma frequency as a measure of the charge density in the system via 

dw Re (Tn{oj) 

In the superfluid phase, this definition excludes the delta function at cj = 0. In the case 
of the superconducting phase transition, where the spectral weight is transferred into the 
delta function at w = 0, the degrees of freedom can rearrange themselves but they cannot 
be lost. The Ferrell-Glover-Tinkham sum rule can also be expressed in the form 

TT 

Ps = — lim [a;Im(Ts(a;)] = / dw [Recr,i(a;) — Recrs(a;)]. (3.17) 

2 a;-5-0 Jq+ 

Here an denotes the optical conductivity in the normal phase, i.e. for some T > as 
the conductivity for some temperature below Tc, and ps is the superfluid density at that 
temperature. The contribution from w = 0 has been separated out explicitly giving rise to 
the term ps determined by ag- According to (3.17), the superfluid density is equal to the 
missing spectral weight, i.e. the difference in the area under the conductivity curve in the 
normal and in the superconducting state, c.f. Figure 13. Note that (3.17) assumes already 
that the translational symmetry is broken, i.e. a 6{oj) contribution or diamagnetic pole in 
the normal phase is absent. It is convenient to define^^ 



(3.16) 


1 

T{ujc) = — / dw [Recj,i(a;) — Re(Ts(a;)] , (3.18) 

Ps Jo+ 

in order to apply the sum rule to the numerically calculated conductivities. Here Uc is 
a cutoff frequency and the sum rule is satisfied if X(oo) = 1. In Figure 14, I{oJc) is 
plotted in the condensed phase for two different temperatures ^'/fi = 0.1, 0.05. As expected, 
liujc) approaches unity for large enough cutoff frequencies. This confirms the sum rule 

more elegant way of regularizing (3.15) and (3.16) is nicely described in [55]: Instead of working with 
the Green’s functions obtained naively from a holographic calculation, which typically do not vanish in the 
upper half frequency plane and on the real axis for |a;| —>■ oo, defining a subtracted Green’s function with 
these problematic contributions removed ensures that the sum rules are valid. These local subtractions 
correspond to the addition of local finite counterterms to the holographically renormalized partition func¬ 
tion. In this way a particular preferred renormalization scheme can be chosen without invoking additional 
requirements such as supersymmetry [56]. An example is described around eq. 12 of [39[, where a local 
ioj/2 term (also present in (3.9)) was removed by a finite counterterm f d*Xy/—^Fij. We thank Martin 
Ammon for pointing us to the latter reference. Note however that whatever renormalization scheme is 
chosen, UV counterterms can only affect the ultralocal terms in the Greens function and hence the UV 
asymptotics of the conductivity. The physical part of the conductivity which we are interested in, and 
which after renormalization should fulfill the Kramers-Kronig relations (3.12) (i.e. causality), comes from 
the current-current two point function at different points in space-time and hence cannot be affected by 
this choice, but must be renormalization scheme invariant. 

^^Additionally, the definition of ps in (3.13) includes the factor of arising from the integration in the 
Kramers-Kronig relations, i.e. pa = 


28 



Figure 13. Visualization of the FGT sum 
rule as explained in the text. The blue 
area indicates the spectral weight which is 
transferred into the zero mode. Note the 
tiny regular contribution which resembles 
a key property of high-temperature super¬ 
conductivity and might be responsible for 
the small possible offset in the computa¬ 
tion of ps, in particular it might account 
for the missing charge density in the super¬ 
fluid phase, i.e. the offset between Us and 
Ps displayed in Figure 10. 



FGT sum rule for g = 6 & k = 0 



FGT sum rule for q = 6 k k — '‘-/V 2 



Uc 


Figure 14. The Ferrell-Glover-Tinkham sum rule holds in the condensed phase for different values 
of ^/^l color coded as in Figure 5. The left panel displays the FGT sum rule in the case k = 0 at 
the temperature '^/^l = 0.1 and Ym (Va*) = 0.3 (0.4,1), 1.5 (0.2,1.4), 3.3 (0.2,1), 6(2.4), where the 
numbers in brackets denote the corresponding values of p/m- The right panel shows the FGT sum 
rule for k = VV2 , Ym = 0.05 and Ym (Ym) = 0.3 (0.3,1 ), 1.2 (0.2), 3 (0.2,1,1.4), 6 (1.4). The integral 
T{ujc) defined in (3.18) measures the missing spectral weight up to the cutoff Wc and is normalized 
so that, if the sum rule is satisfied, X{po) = 1. As can be seen e.g. from Figure 6 for Ym > 8 the 
optical conductivity enters the conformal regime, i.e. for d = 3 -I- 1, cr(w) ^ w, irrespective of the 
existence of the “superconducting gap”. In the conformal regime the normal phase and condensed 
phase optical conductivity becomes identical and thus will not contribute to X{uj). Note that the 
thick black line in the left panel represents the translationally invariant case, p/m = Ym = 
as expected the FGT sum rule fails spectacularly, owing to the coexistence of a ideal metal and 
a superconductor. Thus, in this case the diamagnetic pole in the imaginary part of the optical 
conductivity includes not only the missing spectral weight, but also the ideal metal contribution. 


for the system under consideration and can be seen as a powerful consistency check of the 
holographic model and of the calculation including the numerics. Physically, it shows that 
the charged degrees of freedom of the system are conserved. In particular, it uncovers 
the main obstacle in defining a proper superfluid density in the translational invariant 
system, since the FGT sum-rule as defined in (3.18) does not hold due to the coexistence 
of the normal state ideal metal contributing to the diamagnetic pole, c.f. the black line 
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in Figure 14. Once we turn on our helical structure the “spurious” contribution due to 
momentum conservation are removed from the diamagnetic pole and the FGT sum rule 
confirms the conservation of charged degrees of freedom. 

3.4 Checking Homes’ and Uemura’s relations 

There are two very intriguing relations that were found experimentally, namely Homes’ 
relation [33, 34, 57] and Uemura’s relation [58]. The former is given by 

Ps{T = {)) = Cot)c{Tc)Tc- (3.19) 


whereas the latter reads 


Ps = BTc, (3.20) 

with B being a proportionality constant of units mass in natural units. Uemura’s relation 
is found to hold for underdoped cuprates only, while, as demonstrated in [33, 34], Homes’ 
relation holds for a much broader class of materials. Concerning the units of Homes’ con¬ 
stant, as defined in Section 3.2, ps is given in units of p? and cjdc well as T in units of 
p. Thus, Homes’ constant given by p«/{(jj:icTc) is dimensionless in our unit system. 


Checking both Homes’ and Uemura’s relations by plotting ps against Tc and ctdc^o we con¬ 
clude that Uemura’s relation does not hold in our helical superconducting system. Homes’ 
linear scaling relation, on the other hand, is clearly visible in Figure 15, where we show a 
log-log plot of Ps vs (Jdc^c for various with «: = 0. In the range of 

-Ri4.5, ...,6 and -?al,...,2. (3.21) 

p p 

the relation is linear and extracting Homes’ constant, we find it to be 

Cw 6.2 ±0.3, (3.22) 

Here the uncertainty is not statistical, but refers to the ±5% band bounded by the dashed 
lines in the inset in Figure 15. Intriguingly, comparing our Homes’ constant with the 
experimentally found values [33, 34], after correcting for the factor 8 in our definition of 
the superfluid density, c.f. footnote 14, the helical system seem to be interpolating between 
the dirty limit BCS superconductors with C = ps 8.1 and the in-plane high Tc cuprates 
result C = 35/8 ~ 4.4 [34]. The error bound may be retrieved from [33] by converting 
from the dimensionful constant in units of to our dimensionless unit system^^ 

and yields C = 4.4 ± 0.9. In fact the value of C ~ 6.2 presented in Figure 15, seems to 

the data analyzed in [33] the unit of Homes’ constant is given by 

[rl = [Pd ^ ^ cm~^ 

kDc][T] n-icm-iR fl-iR' 

To convert to a dimensionless unit system used in our holographic system one needs to introduce the natural 
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Homes’ relation for g = 6 & k = 0 



Figure 15. A log-log-plot of the superfluid density Ps with respect to (Tdc(Tc)7c- The color coding 
for Ym is identical to the phase-diagram plots presented in Figure 5, Ya* = ' 4.8, 5.1, 5.4, 5.7 

, whereas values of different p/m are not resolved, except for the outliers where the value of p/^^ is 
explicitly attached to the point. If Homes’ relation holds, the points should roughly fall on a line 
with a slope of unity, according to log(ps) = logC -I- log(CTDcT'c) denoted by the black line. The 
inset shows the value of Homes’ constant C for Ym = 4.5, 4.8, 5.1, 5.4, 5.7. The relation is not 
expected to hold in the limits of p/m —)■ 0 and p//x —)■ oo. There the constant approaches zero due to 
the absence of momentum relaxation and the corresponding divergence of ctdc- These data points 
may be faithfully discarded. Doing so, we see that, in the reasonably applicable range of p/^^ G [1, 2] 
Homes’ relation seems to hold within the dashed lines given by C ~ 6.17 ± 0.31. This value for the 
constant is extracted from a least-squares fit represented by the thick black line in the main figure. 


be almost the arithmetic mean of the two experimentally determined values. Additionally, 
one may compare to the most recent results found for organic superconductors in [57], i.e. 
C = (110 ± 60) again in dimensionful units. Converting to our dimensionless 

Homes’ constant and including the additional factor of 8, we hnd C = 4 ±2.1, which is very 
close to the original result in [33]. 


constants, e.g. for the conversion of the temperature we have 

T[K] = ^.100T [cm-'] , 

which amounts to IK = 0.695cm-'. Similarly, Iff-'cm-' = 4.935cm-' and our final conversion factor 
reads Iff-'K = 3.42983 cm-'. Thus, the values given in [33] are converted by 


(120 ±25) 


cm 

fl-iK 


120 ± 25 
3.43 


35 ±7.3. 


Taking into account the correction factor for our different definition of pa we arrive at C = 35±7.3/g 
4.4 ±0.9. 
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Homes’ relation appears to hold for high values of Y/j > 4.5. This is the regime where the 
Drude peak in the optical conductivity becomes very broad such that the Drude regime, the 
intermediary frequency regime and the conformal regime are shrinking together, indicating 
a mixing of IR and UV degrees of freedom. One observation from the finite temperature 
thermodynamic phase diagram, c.f. Figure 5, is that for such high values of Ym we may hit 
a quantum critical point at a critical value of pY- On the other hand. Homes’ relation, as 
shown in the inset of Figure 15, seems to work over a finite range of p/im beyond the possible 
quantum critical point at T = 0, which clearly calls for further investigation of the zero 
temperature system. 

Alternatively, according to the single scaling argument given in [59] Homes’ relation seems to 
require two competing timescales. In our system the helical lattice introduces an additional 
timescale for momentum relaxation, controlled by Ym; which is very different from the 
diffusive timescale in the original holographic s-wave superconductor c.f. [36], at least for 
small values of Yp- It is compelling to speculate that in the large Yp regime, where the 
applicability of the Drude model may be problematic, these two timescales may become 
almost identical. Let us stress that for a complete understanding of the aforementioned 
scaling relations it is imperative to understand the zero temperature phases of the helical 
system. Nonetheless, the optical conductivity with its broad Drude peak resembles the 
dirty limit BCS superconductors, where Homes’ relation follows naturally from the missing 
spectral weight argument: due to the broad peak, we may think of the missing spectral 
weight area roughly as a square spanned by ctdc and the width of the gap, which is set by 
the universal gap equation at T = 0 to be a number times Tc, see also [60]. We comment 
on this further in our discussion in section 5. 

4 Zero Temperature Solutions and Holographic RG Flows 

In order to solve the system at zero temperature and to understand its zero temperature 
phase diagram and quantum phase transition structure, it is necessary to identify the correct 
infrared geometries. Classifying all possible IR geometries is in general a complicated task 
which can only be done by restricting to certain Ansatze and symmetry requirements, but 
within that class may yield interesting physical insights [61-63]. As a possible candidate 
for such a solution, we now generalize the insulating geometry of [26] to the case of an 
additional massless charged scalar. We want to emphasize that this insulating geometry is 
different from the usual gapped AdS-Soliton geometry: Instead of the holographic direction 
ending at some particular point, the IR of this solution is an anisotropic hyper-scaling 
violating Lifshitz throat. This anisotropy forces the system to be a smectic material, i.e. 
an insulator in the direction of the helix (the x direction), and a metal in the other two 
orthogonal directions. 

For reasons explained in Section 2.1 we work with a non-vanishing Chern-Simons coupling 
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n = (but still m = nip = 0).^^ The solution can be written as a power series in 
with the leading terms being 

re = rco + -\ -, P = Po + pir‘^^^ H-, a = H-, 

18 

e’^i = e’^io r"'/'+ • • • , e’'^ = e^^o + ... ^ u = —r^ + ■ ■ ■ (4.1) 

5 

e^3 = e^3o ^1/3 + ... _ 


The coefficients of this expansion can be expressed in terms of the parameters po, V 20 and 
^ 30 -^^ In particular, it follows from the equations of motion that 


«o — 


Pi = 


36 k sgn(p) e^’^^o ‘^I’so 

5 (Gk^ + q'^Pq — 4 ) ’ 

e^^2o-4t;3o 

( 6 k 2 + q^pl - 4 ) 


2 ’ 


u;o = ^e2^30-^20 ^ 


e^-io 


\/3 (gVo - 4) e3^2o-2t;3o 

T(6k 2 + g2p2 -4^ - ' 


1 

2 


\p\ e^30-«20 ^ 


(4.2) 


This fixed point describes a cohesive IR geometry with a superconducting order parameter 
turned on. Note that the charged scalar is not subleading or leading compared to the 
original geometry without it, it rather has the same IR behavior as the helix field w. 


In order to understand whether this fixed point is stable under perturbations, we follow 
[26] and calculate power law perturbations around (4.1) by writing 

w = Wo + + CwT^), p = po + pir'^^^il + Cpr^), a = aor^^‘^{l + Car^), 

18 

Vi = Ulo + log(r“^/3) = .^20 + log(r^/®) + C2r^, U = —r^(l + cur^), (4.3) 

5 

^^3 = '^30 + l0g(r"/3^ 

The equations of motion are linearized in the perturbations and solved to leading order in 
r. All scaling exponents and the corresponding eigenvectors for the radial perturbations 
around (4.3) are listed in Appendix C.3. In summary, we find that the condensed insulating 
solution (4.3) does not show any condensation instabilities in which some of the IR operator 
dimensions <5 violate the Breitenlohner-Freedman bound by becoming complex. Instead, we 
find two IR irrelevant deformations, i.e. deformations with explicitly positive exponent, 
namely the 5+ mode of point 6 and 7 in Appendix C.3, 

61 = ^ (Vl 45 - 5) , 62 = ^ (VM - 5) . ( 4 . 4 ) 

These perturbations will be useful in generating the RG flows to the UV by shooting 
numerically from the IR fixed point perturbed with these deformations (c./. Appendix D 
for more details). These exponents are the same as the ones found in [26]. As will be 

330ne reason is that it seems harder to find IR scaling geometries for non vanishing masses. We plan to 
return to this question in the near future [37]. 

3^There is an additional free parameter, namely the expansion point r+, which has been set to zero for 
simplicity. It can be reinstated by shifting r —>■ r — r+. 
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explained in detail in Appendix D, these two modes are sufficient to generate the two- 
parameter family of zero temperature RG flows labeled by the chemical potential /^/p and 
the lattice strength ^/p. 


Besides the above superconducting IR geometry our model admits, at least for large enough 
Chern-Simons couplings such as our choice k = two other IR fixed points: For a 

vanishing charged scalar, there is a metallic AdS 2 x fixed point dominating for larger p 

[26], whose geometry including perturbations reads 

w{r) = Cy,r\ p{r) = Cpr\ 

vi{r) = i;io(l Cl/), V2(r) = i;2o(l + C 2 /), 

V 3{ r ) = 1^30(1 + C3/). 

The metallic AdS 2 x geometry has several deformation exponents, which are spelled 
out together with the corresponding eigenvectors in Appendix C.l. Here we focus solely on 
the condensation instabilities. In particular, there are two modes corresponding to scalar 
condensation in this AdS 2 near horizon geometry, with scaling exponents 

S± = ^ (-3 ± y^3(3-2g2 + m^2)^ , (4,6) 

If these exponents become complex, the charged scalar destabilizes the AdS 2 , and the 
system presumably flows to the superfluid IR geometry above. Furthermore, there are two 
exponents connected to the condensation of the helix field. 


a(r) = 2\/6r(l -h CaV^), 

C/(r) = 12/(1 +Cl//), (4.5) 




1 ^ ((m2 + 3) + p2 e-21^10 _ 2\/6Kpe'"^o-2v2o ^ 

2 ^ 2/1 


(4.7) 


Because the lattice is explicitly introduced, the crucial aspect for the physics is now whether 
the exponent becomes relevant. This happens when <5+ < 0, see [26]. In that case the 
system will flow to the insulating geometry. If the exponents become complex, then the 
AdS 2 geometry can spontaneously destabilize to the insulator, but we will not consider this 
particular case. The insulating geometry of [26], given by (4.1) with the charged scalar p 
switched off, is also unstable towards condensation of the charged scalar within the system 
(2.2), although in a slightly different way. Analysing the radial perturbations for the case 
of vanishing scalar mass rup = 0 one finds an additional mode for the charged scalar alone. 


5p = co + cir 


5/3 


(4.8) 


If the charged scalar had a non vanishing mass, its exponents would change from (4.8) to 

^®For Chern-Simons couplings smaller than the critical value Kc ~ 0.57 another unstable IR scaling fixed 
point appears [26], which complicates the phase structure at zero temperature. Here we discuss only the 
simpler case of large k. 


(4.9) 
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Note that the IR dimension of the charged scalar in the insulating background of [26] is 
independent of its charge q, due to the cohesive nature of the extremal horizon. In the 
regime 

-4<m2<-^, (4.10) 

the charged scalar obviously violates the IR Breitenlohner-Freedman bound while preserving 
the UV Breitenlohner-Freedman bound, and the condensation mechanism will be analogous 
to the metallic AdS 2 x IR^ case. On the other hand, for the massless case (4.8), no conden¬ 
sation instability is found. In this case, condensation can still happen thermodynamically if 
the condensed zero temperature RG flow obtained from the IR geometry (4.1) has a lower 
free energy compared to the uncondensed one (Eq. (4.1) with po = 0)- We numerically con¬ 
structed the holographic RG flow geometries up to the asymptotic AdS boundary for both 
the insulating and superconducting fixed points for a certain range in parameter space, and 
confirmed that they have lower free energy. For completeness, we collect all the operator 
dimensions and the corresponding modes for each fixed point in Appendix G. 

5 Discussion and Outlook 

In this work we analysed the transition to s-wave superconductivity in an anisotropic five¬ 
dimensional holographic model with a helical Bianchi VIIq symmmetry. This corresponds 
to a 3+1 dimensional field theory in the presence of a helical lattice [26]. The advantage 
of this model is that it allows us to cleanly separate the IR dynamics in the system. This 
is hard to identify in the simplest holographic superconductors for two reasons: Due to 
translation invariance there is already in the normal phase a delta peak at zero frequency in 
the conductivity. In the superconducting phase this mixes with the protected fluctuations of 
the order parameter. Secondly, most well-known examples of holographic superconductors 
are accompanied by a remaining gapless Lifshitz sector in the IR that mixes dynamically 
with the order parameter physics. This is especially so at finite temperature. We improved 
on the former point by explicitly breaking translation invariance along one of the field theory 
directions using the above-mentioned Bianchi VIIq helix, and on the latter by using the fact 
that this model (2.2) has an anisotropic insulating ground state [26]. We established that 
this model indeed undergoes a superconducting transition at low temperatures. Studying 
the optical conductivity we can see that the IR dynamics is more cleanly controlled by 
order parameter physics. This allowed us to extrapolate to a first holographic example 
where Homes’ relation holds. Let us discuss the physics of each of these points. 

5.1 Phases at Finite and Zero Temperature 

The phase diagram of the holographic helical Bianchi VIIq lattice model is quite rich and 
this is reflected in the ways it approaches superconductivity. For large enough charge q 
of the scalar order parameter, both the insulating phase at small helix pitch as well as 
the metallic phase at large helix pitch are unstable towards condensation of the charged 
scalar. The second order mean field superfluid transition typically happens at a critical 
temperature rc(Y/i, p//x), but the data in Figure 5 suggests that a quantum phase transition 
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between condensed and uncondensed phases is possible for larger values of ^/^l, similar to 
the situation in a recently investigated axion-based system [30]. 

The curious aspect is that the critical temperature does not have a monotonic behavior 
as a function of the helix parameters. Naively the presence of a lattice should form an 
obstacle for s-wave superconductivity. This is true at very small helix parameters. There 
Tc decreases compared to the translationally invariant system. However, for a given value 
of the amplitude A there is a critical value of the helix pitch p beyond which Tc starts to 
rise again. In the presence of the Chern-Simons coupling, n = VV 2 , Tc can even increase 
beyond its isotropic value for very large p/^l. In the absence of the Chern-Simons coupling, 
arguably, the tendency to return to its original homogeneous and isotropic value for large 
v/^i can be understood as the effect of the helix diminishing if it rotates too fast around 
the x-axis: The valleys between the maxima become so narrow that they do not influence 
the condensation dynamics any longer, and homogeneity is approximately restored. It is 
an open question whether all observables return to their homogeneous values at large p/^l, 
and at which rate. 

In neither case, however, is the physics behind this behavior of Tc very clear. There is a 
strong indication, on the other hand, that it is correlated with the zero temperature ground 
state of the system in the normal phase. We did not construct all of these, but one can 
infer from the hnite temperature optical conductivity qualitatively whether the true ground 
state is insulating or conducting, see Figure 16.^® These do indicate a second insulating 
phase occurs at large helix pitch p or equivalently small helix wavelengths. One now sees 
that there is a rough correlation between high Tc with an anisotropic insulating ground 
state in the normal phase and low Tc and a metallic ground state in the normal phase. The 
correlation is not exact, however. Clearly, an independent analysis from thermodynamic 
quantities as well as a complete calculation of the zero temperature phase diagram is re¬ 
quired to establish this concretely and unambiguously decide the fate of this new insulating 
phase.The correlation of the behavior of Tc with the zero-temperature normal phase 
ground states indicates that the naive insight that homogeneity is approximately restored 
is probably incorrect, as then the system is expected to be in a conducting rather than an 
insulating phase. A brief investigation into the possible zero temperature ground states, 
allowed us to construct an IR geometry dual to the superconducting phase based on the 
original insulating solution of [26], c.f. Figure 1. Interestingly, the charged scalar shows the 
same approach to the IR as the helix held, indicating that they might be able to compete 
in quantum phase transitions. We analysed the static radial perturbations around these 
three IR hxed points, in order to understand which RG hows between them are allowed. 
The situation is summarised in Figure 1: The metallic AdS 2 x IR geometry behaves 
conventionally. It can be unstable towards either the insulating state and/or superconduc¬ 
tivity [41]. At the same time the condensed superconducting IR geometry we constructed is 

^®Note that the finite temperature solution is uniquely determined from the boundary conditions. It 
therefore already knows whether it originates from an insulating or a conducting zero-temperature geometry. 

^^We thank Aristomenis Donos for discussions on this point. 
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Figure 16. The nature of the zero-temperature ground state from the finite temperature conduc¬ 
tivity. The surprise is that for fixed helix amplitude A the system transitions from an insulating to 
a metallic and then back to an insulating phase. For p/^^ ^ 1 and Ym ^ li we expect a metallic 
phase designated by the shaded green area due to the fact that momentum relaxation is removed 
in the limit where either of these parameters vanishes. This part of the metallic phase could not 
be distinguished from the insulating phase since we are using very coarse measure to determine 
the nature of the ground state i.e. a qualitative measure of the conductivity. The thick blue line 
denotes the location of the minimal critical temperature T/”'" extracted from Figure 5. In the case 
At = 0, shown in the left panel, this minimum tracks qualitatively the metallic phase inferred from 
the conductivity, whereas for n = i/Va the minimal critical temperature T/”’” is invariant under 
changes in p/m s-nd Ym- Note that at high values for Ym th® critical temperature is very low and 
thus our numerical code cannot reach anymore. Homes’ relation holds in the region marked 
by the white dashed box. 



nicely stable, indicating that it is the true ground state [26]. The insulating IR geometry is 
indeed unstable towards superconductivity, but curiously not for the mass of the scalar field 
considered here. We suspect, however, that in this case the superconducting IR geometry, 
is still the thermodynamically preferred ground state, i.e. the state of lowest free energy. 
The insulating but not superconducting geometry of [26] is hence dynamically stable, but 
thermodynamically unstable. This would indicate that they are separated by a first order 
transition. We will support this claim by an analysis of the thermodynamics and transport 
at zero temperature in a forthcoming work [37]. 

5.2 Transport 

In our system, the linear momentum relaxation introduced by the Bianchi VIIq structure of 
the geometry allows us to reliably analyse the physics behind the low-frequency transport 
properties of our system. Our computation reveals that the superconducting system is well 
described by a two-fluid Drude model at small frequencies in the regime of weak momentum 
relaxation Ym <C 1, a fact also observed in the models of [29, 30]. On the other hand, in 
the regime of stronger momentum relaxation, Ym ~ Ij the two-fluid Drude model seems to 
work less and less well, again similar to [29, 30]. 
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It is the absence of a diamagnetic pole in the normal phase, that allows us to reliably extract 
the superfluid density from the i/w pole in the imaginary part of the optical conductivity 
in the superconducting phase. Naively the cohesiveness of the superconducting phases in 
our system at zero temperature literally forces all charge carriers in the system to condense 
into the matter fields outside the extremal black hole horizon.^® Since the U{1) gauge 
field sourcing the helix itself is not charged under the ‘charge’ U{i) and since the Chern- 
Simons term in (2.2) only induces currents but not a charge density for the ‘helix’ 17(1), one 
would expect all the charge density (as carried by the ‘charge’ 17(1)) in the system at zero 
temperature to be carried completely by the charged scalar dual to the superconducting 
order parameter. The charge simply has no other place to go within this system. By 
comparing the charge density in the superfluid phase at very low temperatures to the 
superfluid density as extracted from the optical conductivity (c./. Figure 10), we showed 
numerically that these indeed become identical for small lattice strengths ^/^l 1. In 

more detail, we found that when the Drude model is a good approximation to the optical 
conductivity at small frequencies, i.e. when Y/i <C 1, the difference between the charge 
density at small temperatures and the superfluid density is accounted for by the additional 
spectral weight residing in a small residual Drude peak at low frequencies. The fact that we 
are analyzing the system correctly is confirmed by the fact that the Ferrell-Glover-Tinkham 
sum rule holds (c./. Figure 14) when crossing the phase transition from the normal to the 
superconducting phase, i.e. that there is no missing spectral weight in our system. This is 
to be contrasted with the translationally symmetric case (the black line in the left panel 
of Figure 14), in which the sum rule fails spectacularly due to the non-accounted spectral 
weight residing in the 5{uj) poles related to momentum conservation in both the normal 
and superconducting phases. 

For weak lattice strengths Ym Ij the model considered here is therefore in several respects 
under better control compared to the simple holographic superconductor analyzed e.g. in 
[36]. A puzzle appears for stronger momentum relaxation Ym ~ 1- Now the charge density 
and superfluid density do not approach each other up to the small temperatures probed in 
our numerics, and the difference cannot be accounted for any more by a normal fluid Drude 
component. Even so, the FGT sum rule continues to hold. A possibility is that some of 
the low frequency spectral weight gets transferred to intermediate frequencies rather than 
the superfluid pole. Another explanation for the failure could be that our identification of 
Ps from the optical conductivity does not match with the superfluid density as calculated 
from the transverse response via the magnetic/London penetration depth, which is another 
important cross-check [37]. This deserves further study. Based on the new insulating phase 
found at high p/mj ^ distinct possibility is that the phase structure of the system is more 
complicated and interesting for strong lattice potentials at high Yp- It is precisely in this 
novel regime where we can find a region in parameter space in which Homes’ relation is 
valid to a good accuracy. 

■^^Note that a translationally invariant cohesive phase still can have a zero frequency delta function due 
to the presence of the charge density, as explicitly shown in [64]. 
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5.3 About Homes’ Relation 


Homes’ relation [33, 34] - the experimental result that for high Tc superconductors as well 
as conventional BCS superconductors, there is a universal relation of the form 

PsiT = 0) = CauciTc)-Tc, (5.1) 

with a nearly universal, material independent constant C — would follow naturally from 
an argument based on the shortest possible time scale in a strongly coupled quantum crit¬ 
ical state [35, 51], a so-called Planckian dissipator. The basic idea is that in a quantum 
critical system there is no other scale than the temperature, and hence the relaxation time 
responsible for the finite electric conductivity in the system must depend on the inverse of 
the temperature alone, up to numerical factors of 0(1). Holography is unique in its ability 
to model interacting quantum critical systems and the notion of Planckian dissipation is 
clearly visible in the universal holographic result for the ratio of shear viscosity over entropy 
density [65]. A holographic foundation for Homes’ relation therefore has the potential to 
confirm that a similar universal mechanism is at work in superconductivity.^® 


In our system, however, the timescale is not the intrinsic timescale associated with the 
strongly coupled dynamics of the system. Instead it is the external momentum relaxing 
timescale introduced by the lattice. This is evident from the validity of the Drude re¬ 
sponse, where ctdc ~ 'Tmomentum- The relaxation time scale that controls the low-frequency 
conductivity is thus a function of the lattice strength and helix pitch. 




T 


(5.2) 


We have extracted an explicit factor of T such that the function / is dimensionless. A 
Planckian dissipator would have a mostly constant function / of order 0(1) in a CFT. 
Clearly, for weak momentum relaxation, the relaxation time is typically very large, tT ^ 1, 
non-universal and far from a Placnkian dissipator. 


The breaking of translation invariance and the introduction of an external scale is, however, 
important in studying Homes’ relation with gauge/gravity duality. As we emphasized, 
in a translationally invariant system there is already an infinite (5-function contribution 
to the DC-conductivity. To extract the superfluid density reliably, one needs to resolve 
this either into a Drude peak behaviour for weak breaking or to something beyond [67- 
69]. Additionally the system should be in a cohesive phase at low temperatures in order 
to prevent the existence of additional IR charged degrees of freedom, hidden behind the 
extremal horizon. It was found in [1, 4, 52] that a modulated chemical potential is not 
sufficient to realise Homes’ relation in the simple model considered there. In view of the 
above, these works probably did not access the regime of strong momentum relaxation 
that could potentially make Homes’ relation work. For example, in [1, 4, 52] the two-fluid 

course, though actual transport in high Tc cuprates does have Planckian dissipative features, it 
cannot be a pure quantum critical state, see e.g. [59, 66]. 
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Drude model works well for all lattices considered, pointing to a regime in which momentum 
relaxation is weak. Instead, in this work a clearer picture emerges: Due to the cohesive 
nature of the superconducting ground state, as well as the broken translation invariance in 
the helix direction, we were able to eliminate, respectively, the additional charged degrees 
of freedoms in the IR and the zero frequency delta function, hence obtaining a clearer 
account of the reshuffling of charge density as well as spectral weight during the superfluid 
transition. As noted already in [36], to successfully analyse Homes’ relation in a holographic 
model it is essential to be able to keep track of all charged degrees of freedom at low/zero 
temperatures. Since our helical lattice model admits an insulator/superconductor transition 
with a cohesive phase at low temperatures, this setup allows for a fresh look at Homes’ 
relation, addressing both of these points. 

Figure 10, which shows the agreement between the superfluid density and the total charge 
density in the limit T —>■ 0, illustrates this cleanly, together with our verification that 
the Ferrell-Glover-Tinkham sum rule holds. The validity of the sum rule is of particular 
importance, as it can be used to obtain Homes’ relation if the underlying system is a 
Planckian dissipator [36]. With these results as well as (5.2), (5.1) can be reformulated as^^ 

= C'(A/;.,p/^)/(T,//r, ...). (5.3) 

nn{T = Tc) 

with C the coefficient in Homes’ relation. In simple systems, with a single species of charge 
carriers and a reasonable gap in the superconducting phase, the charge density at zero 
temperature in the superfluid phase is approximately equal to the charge density at the 
critical temperature in the normal phase. From (5.3), Homes’ constant must hence vary 
inversely with the function / parametrising the momentum relaxation time scale, 

C ~ 1//. (5.4) 

In a bona fide Planckian dissipator, / would be a universal constant of 0(1), and is in this 
way seen to explain Homes’ relation. We also investigated the LHS of (5.3). Figure 17 shows 
that numerically the charge density in the system does not vary much between T <C Tc and 
Hence in combination with our results from Section 3.2 that the charge density and 
superfluid density are approximately equal at low T, also the ratio of charge densities in 
(5.3) is seen to be close to unity. This is so in the weak momentum relaxation regime 
where A/^ <c 1, but also for rather large i-e. strong momentum relaxation in which the 
Drude model is not readily applicable, in this case Homes’ relation may hold. As we found 

^°Note that compared to the standard condensed matter notation we absorbed the effective masses and 
nnmerical factors in the respective phases (which are not directly accessible in holographic models) into n 
by n„e^/m* i—>■ n„ and rise^/m* i—>■ Us, i.e. we used ctdc = nr and ps = Us- This redefinition gives n„ and 
Us units of (length)”^, which coincides with charge and not number densities in five dimensions. 

^^Figure 17 does indicates that for most choices of Ya* and p/m our system has slightly (~ 5 ... 10%) more 
charged degrees of freedom in the superfiuid than in the normal phase. While in the normal phase the 
system is anisotropic, the system becomes more isotropic in the superfiuid phase due to the isotropic s-wave 
order parameter. An explanation for the rise could be that charged degrees of freedom from the metallic 
transverse directions also contribute to the charge density in the helix direction. 
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in Section 3.4, explicitly computing the functional relation between ps and <tdc(2c)Tc we 
indeed find that in our holographic superconductor in a helical lattice Homes’ relation holds 
with coefficient 


Cw 6.2 ±0.3 


(5.5) 


in the range of parameters 

— « 4.5,..., 6, and 2, for k = 0. (5.6) 

p p 

It can, however, not be explained as a consequence of the simple ratio of charge densities 
at Tc and T = 0. This is also the regime where a difference between ps and n* opens {c.f. 
Figure 10), and the simple relation between (5.3) and (5.1) breaks down. 


As we now argue, the fact that Homes’ relation does hold in the parameter range given above 
may on the contrary be due to a strong deviation from the simple Planckian dissipation 
behavior.We may parametrize the difference between ps and Ug as 


Ps{T = 0) = ns{T = 0) ± 5ps{T = 0). 


(5.7) 


Furthermore, let us assume^^ the DC conductivity is still roughly proportional to the charge 
density times a relaxation time, 

o-DC = nnT. (5.8) 

Substituting Homes’ relation on the LHS of (5.7) together with (5.8), and using that the 
ratio (5.3) is close to unity, we arrive at 


c/(r„A/^,p/^,...) 


1 + — 0 , . . . ) 


(5.9) 


Then, for Homes’ relation (5.1) to be valid and C to be a univeral constant of 0(1), we see 
that the gap 5ps and the function / cannot vary independently from each other, but must 
conspire. As 6ps varies with the parameters, / cannot be constant, and hence the system is 
not a Planckian dissipator. A more detailed investagation of the function / directly as well 
as other relaxational scales is hence of great interest, as well as how the gap 5ps behaves 
in the regime of strong(er) scattering. Since this is no longer Drude physics, it requires a 
more detailed response analysis to extract these. We do note that the regime of p/^l and 
A/^ where we find Homes’ relation to hold is near the apparent insulator-metallic quantum 
phase transition in the normal phase, see Figure 16. It will hence be very interesting to 
investigate in more detail the behavior of the above quantities in the zero temperature 
ground states [37]. 

are grateful to E. Kiritsis for pointing out this possibility. 

^®We make this assumption in hindsight of the parameter regime in which Homes’ relation is valid (see 
below) to be not too far into the regime of strong scattering, such that this Drude-like approximation still 
should yield reasonable results. Also, it is reasonable to assume that a system conducts better if more 
charged degrees of freedom are present. 
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To conclude, we would like to emphasize that in our model we have found an example of 
the validity of Homes’ relation in a holographic model with strong momentum relaxation. 
This is similar to the experimental result [33, 34] in which the c-axis high cuprates as 
well as the dirty limit BCS superconductors, both materials with strong scattering and 
momentum relaxation, follow the same Homes’ relation. As explained above, in the strong 
momentum relaxation regime there is a nontrivial difference between ps and Ug at low 
temperatures, whose origin needs to be investigated in more detail [37]. An important 
consistency check will be to compute the superfluid strength not from the longitudinal 
response, i.e. the plasma frequency and the Drude weight, but from the transverse response 
by determining the magnetic/London penetration depth. This requires to solve for the 
transverse propagator by finding a solution at small but non-zero momentum. Indeed, an 
analysis of the dynamical conductivity axx{^, k) and possible finite momentum instabilities 
[70] will allow us to determine the spectrum of quasi-particle excitations, and also whether 
the superconducting phase (4.1) is really the thermodynamically preferred ground state in 
our system. There is another interesting aspect to consider, that is the case of spontaneously 
generated helical ground states. Similar to the case of spontaneous generated charge density 
waves [14, 16, 18, 20] this might even dynamically fix the preferred helix pitch p, leaving 
fewer free UV parameters. If Homes’ relation hold for all such spontaneous helix models, 
then one would have a more satisfying reason to explain its universality. At this stage, all 
these results are still preliminary and for k = 0 only, and we plan to analyze the question 
under which exact conditions Homes’ relation holds and related questions, in a future work 
[37]. 
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Figure 17. The charge density as a function of for p/ti = 0.4 and for various values of Ya*) 
using the same color coding as in Figure 5: 0.3, 0.6, 0.9, 1.2, 1.5, 1.8, 2.1, 2.4, 2.7, 3, 3.3, 3.6, 
3.9, 4.2, 4.5, 4.8, 5.1, 5.4, 5.7, 6 in the case of /t = 0 (left panel) and 0.3, 0.6, 0.9, 1.2, 1.5, 1.8, 
2.1, 2.4, 2.7, 3 for k = i /%/2 (right panel). Curiously, the charge density at the phase transition 
seems to be minimal and hence ns(T < T^) > n„(Tt,). Note that this is more prominent at vanishing 
Chern-Simons coupling k. The decrease of the charge density in the normal phase can be attributed 
to the effects of the insulating phase, i.e. there are less charged degrees of freedom available for 
lower temperature. However, in the superconducting phase the system turns completely into the 
condensed state which is growing with lower temperatures. Therefore, the superconducting degrees 
of freedom are increasing for lower temperatures and may be drawn from the metallic directions, 
due to partial restoration of the isotropy. 
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A Equations of Motion for S-Wave Superconductors on a Helical Lattice 


The equations of motion for the metric fields U, vi,V 2 ,V 3 , the Maxwell fields A and B, and 
the scalar field p following from the action (2.2) are 


0 = a" + a' {y'l + + Atpe 
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(A.l) 

The third equation is first order and originates from the rr-component of the Einstein 
equations. There is another second order equation 


1 


+ 4= + i e-"“ + V? + v? + v'i + 4 + 4 + 4' = 0, 


(A.2) 


which follows from the above equations. Therefore, the first order equation is a constraint. 
In the above equations, we have made use of the C/(l) symmetry associated with to 
choose p real. 


A.l Linear Response for S-Wave Superconductors on a Helical Lattice 

We consider linearized fluctuations around the background solution. Therefore, we write 

9fiu = 9p,u + ^^lv, A^ = Mp + Ajj^, iA ll 

= + P = P^ + P^, 

where the fields with a superscript ‘6’ denote the background solutions, h denotes the metric 
fluctuation, and fields with a superscript ‘/’ denote the matter fluctuations. The metric 
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and the vector fields are expanded in the basis (dt,dr, wi,ti; 2 , W 3 ). The background fields 
are r-dependent only and are written in the Ansatz of Eqs. (2.13), (2.9) and (2.10). The 
fluctuation fields are chosen to depend on r and t because the retarded Green function 
leading to the conductivity is evaluated at zero spatial momentum, Eq. (3.1). Expanding 
the action of Eq. (2.2) to second order in the fluctuations, we obtain an action Sq that 
determines the linearized equations of motion. The term linear in the fluctuations vanishes 
upon use of the equations of motions for the background fields. Analyzing the action Sq, we 
can determine which fields couple to each other. The result of this analysis is summarized 
in Table 1 on page 47. The block of fluctuations which contains and decouples from all 
other fluctuations is 

(^A = A{,B = Bl,hti,h23,hri'^ . (A.4) 

It is more convenient to work with the fields 

£ = (/)^^ hti = hti, (A.5) 

-T = h23 = h23, (A. 6 ) 


instead of hti and / 123 . The reason is that £ and T" have a finite limit for r —>■ 00 , whereas 
hti and /123 are proportional to for large values of r. The equations of motion for the 
fluctuations are obtained by varying the action Sq. After variation, we can set hri = 0 
choosing radial gauge in which all radial field components vanish. Carrying out a Fourier 
transform of the time coordinate, i.e. choosing a harmonic time dependence e“^^*, we 
obtain the following linearly coupled ordinary differential equations in r for the fluctuation 
fields A, B, £, and T: 
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The first three equations are second order. The last one - the constraint - is hrst order and 
originates from the metric mode hri after choosing radial gauge. There is a forth second 
order equation, 

0 = T" + £' (3u'i + u' + v'^) + e-2"i a'A' + 

— g-2(^i--y2+i'3) ^ 

+ i^ ^ _ j-e-2^1 + j-e-2(^i-^2+^3) ^ ^ (A.8) 

which follows from the above equations and the equations of motion for the background 
helds. The background equations of motion have been used to eliminate second order 
derivatives of background helds in the above equations. 
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Table 1. Couplings between the fluctuation flelds. Couplings between two flelds in the quadratic action Sq, which determines the linearized 
equations of motion, are indicated by •. The block of coupled fluctuations containing A( is marked in red. 












B Asymptotic Expansions 


Asymptotic expansions of the fluctuation fields are computed near the thermal horizon Vh 
and near the boundary. The expansion around has the leading terms 

A={r- (4 + A>l(r - r,) + • • • ) , 

B = ir- + 4(r - r;,) + • • • ) , 

\ / (B 1) 

£ ={r - + E^ir - ruf • • • ) , 

F={r- (4 + 4(r - r;,) + • • • ) . 

Infalling wave boundary conditions have been chosen, which lead to the retarded Green 
function. The expansion has three free parameters which are chosen to be Ag, Bq, and Fq. 
All remaining expansion coefficient can be expressed in terms of these three parameters by 
means of the equations of motion. We can choose three linearly independent points in the 
{Al 4, 4) -space to define initial conditions for a numerical integration starting from a 
point rh + 5 with 6 numerically small. In this way, we obtain three linearly independent sets 
of solutions. These can be linearly combined to satisfy three conditions at the boundary. 
Two of these conditions are given by requiring that gauge invariant fields are built from B 
and T", and that these have no source at the boundary. The third condition corresponds to 
the normalization of the solution. At the boundary, a double expansion in i/r and log(^)/r is 
carried out. The leading terms of the expansion^^ are 


A = A\^+ ^2 + \og{r)Al/2 ^ 




24 + log(r) [(^2 - p2) Qb ^- 2pF^)] 


2r2 


^ 4-pAlog(r)(pAi^^i^Bg)/4 

A = -^0 H- —A -^ 


(B.2) 


J- = 4 + 


(u;2 - V) 4 4 + log(r) [{co^ - /16 








This expansion has seven free parameters, namely Ag, A2, Bq, B2, Eq, Fq, and F^. The 
coefficient E\ and higher order coefficients can be expressed in terms of these parameters. 


B.l Residual Gauge Transformations and Physical Degrees of Freedom 

In order to determine the physical degrees of freedom corresponding to the fluctuation fields, 
the residual gauge transformations of the radial gauge and their action on the fluctuation 
fields are worked out. The physical fields are those being invariant with respect to a residual 
gauge transformation. To determine them, we follow a similar calculation carried out in [48] 
in the context of the holographic p-wave model. The gauge transformations of the action 

^“^For the sake of clarity, the shift parameter a is set to zero here. It can be reinstated using the 
transformation of Eq. (2.25) and additionally F 4 —> F 4 — 
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of Eq. (2.2) with a massless helix field and vanishing Chern-Simons coupling are diffeomor- 
phisms —7- x'^ = x^ — U{1) transformations of and 17(1) transformations of 

B„ 


j^. Their infinitesimal action on the fluctuation fields is given by 


6Al = + idyAl)j:- + d,A, 

5Bi = + (5.i?;))s"+ d^T, 

5pf = (9.p'’)S" + iqkpK 


(B.3) 


The vector parameterizes the diffeomorphisms, and the scalars A and T the 17(1) trans¬ 
formations. The transformations depend implicitly on the background metric. Since we 
work at zero spatial momentum, we can focus on r- and t-dependent diffeomorphisms and 
U (1) transformations. Furthermore, we work in frequency space assuming a harmonic time 
dependence of A and T. The residual gauge transformations of radial gauge are 

those satisfying 


Shrn = 0 , 6AI = 0 , 6B^ = 0 . 


(B.4) 


Using the transformation rules of Eq. (B.3), a system of differential equations for S, A and 
r is obtained. It has the solution 


Tit = Kt — 

Ty = Ky, 

A = IK'pto 


dr 


1 17(f)3/2’ 

r dfa(f) 

/l [/(f)3/2 ’ 


Tr = KrVU, 

Tz = Kz, 


Ea; — Kx, 


r = ATrr, 


(B.5) 


with Kt^ Kr, Kx, Ky, Kz, K\, and being constants. Having determined the residual 
gauge transformations, we can write down their action on the fluctuation fields: 


5^ = 0, 

6hti = -ie^"^ KxUJ, 


5B = —Kxpw, 

6h23 = Kxp{e^^^ - e-2"^). 


(B.6) 


The gauge transformation of A vanishes and, since e^’^^ = and e^'"® = for large values 
of r, the metric fluctuation h 23 is also gauge invariant at the boundary, where the Green 
function is read out. hti and B can be combined into the gauge invariant field 


G = —iujB -|- wpe 


h 


ti- 


(B.7) 


We therefore have three physical fluctuation fields, namely A, Q, and F = /i 23 - The 

field £ = hti is not gauge invariant and does, therefore, not carry dynamical degrees 

of freedom. 


C Radial Perturbations of Zero-Temperature Fixed Points 

We list below the radial perturbations (IR operator dimensions) and their corresponding 
eigenvectors in the three IR geometries found in our system (2.2). The below notation for 
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the eigenvectors v is connected to the notation in (4.1) respectively by 

(c(7,Ci,C2,C3,Ca,C^,Cp) = CV = c{vi, V 2 , V 3 , V^, V 5 , Vq, V 7 ). (C.l) 

Note that the definitions of the {cjj, ci, C 2 , C 3 , Ca, c^, Cp) slightly change in (4.1), which has 
to be taken into account when using the here-quoted results. In the below tables, M denotes 
an IR marginal, R an IR relevant, and I and IR irrelevant mode. 

C.l Metallic AdS2 x Fixed Point 

Within the system (2.2) including the charged scalar, the insulating geometry of [26] has 
the following radial perturbations: 

1. x-rescalings (M): <5 = 0 and v = (0,1,0,0,0,0, 0), 

2. Combined y/z rescalings (M): 5 = 0 and v = (0, 0,1,1,0,0,0), 

3. Constant shift of the chemical potential (R): 5 = —1 and v = (0,0,0,0,1, 0, 0), 

4. Combined Vi mode (R): 5 = —1 and v = ^0, —2'^, 1,1, 0,0,0^, 

5. Mode in the blackening factor (R): 5 = —1 and v = (1,0,0,0,0,0, 0), 

6. Constant shift of the blackening factor (R): 6 = —2 and v = (1,0,0,0,0,0,0), 

7. Combined gauge field and geometry mode (I): 

5 = 1 and v = (^-^, 1, 0, 0^ , 

8. Scalar condensation mode (if real: 5_ (R), 5+ (I) for > 2q‘^ and (R) for < 

2q‘^): V = (0, 0, 0, 0, 0, 0,1) and 5± = g ^—3 ± ^ 3(3 — 2q‘^ + m^)^. Note that for 
the massless charged scalar used in this work, the 5+ mode is always relevant if the 
exponents are real, i.e. for < 3 / 2 . 

9. Helix condensation mode (if real: 5_ (R), 5+ (I) for — 2\/6 pn > 0 

and (R) for m? — 2\^e~‘^^° pn < 0): 

5± = —^ =F \ + p"^ — 2\/6 e“^io pn )^ and v = (0,0, 0, 0, 0,1, 0). 

Note that for the massless helix field chosen in this work, 5+ in tendency will be 
irrelevant if is large, but the actual value of viq is of course given by UV data. 

10. Additional geometry mode (5_ (R), 5+ (I)): 5± = —5 (^1=F -^1-1- ^ and 

V = (0, 0, —1,1, 0, 0, 0). Note that while for the other modes considered above the 
U"{r) equation enforces C2 = C3, this is not the case for this mode, since here the 
U"{r) equation is proportional to the quadratic polynomial in 5 which vanishes for 
the solutions 5± considered here, and hence is automatically fulfilled to first order in 
the perturbations. 
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Note that the background has V 20 = U 30 in the deep IR. Interestingly, only the last mode 
breaks this as one flows up to the UV. Also, all modes are either marginal, have 5 = —1, or 
come in pairs which sum up to —1. This is obvious for the modes 8,9,10, but in fact modes 6 
and 7 are also a pair arising from the polynomial 6 ^ + 6 —2. Finally, note that several modes 
(5,6,7) contribute to the perturbation of U{r) and hence can contribute to the temperature 
perturbation. This means that one must continue these IR perturbations to the UV by 
constructing the perturbed RG flow in order to understand the different contributions of 
these modes. We will come back to this in future work [37]. 

C.2 Insulating Fixed Point 

By switching off the scalar by setting po = 0 in the zero temperature solution (4.1), (4.2) 
implies pi = 0 and hence the Ansatz (4.3) would have no radial perturbation for the charged 
scalar p. In calculating the above modes we hence used a slightly different Ansatz for p 
compared to (4.3), namely 

P = Po + +Cpr\ (C.2) 

'-V/-' 

=0 

Within the system (2.2) including the charged scalar with mass nip, the insulating geometry 
of [26] {i.e. (4.1) with po = 0) has the following radial perturbations: 

1. Blackening factor mode (R): 6 = —S/s and v = (1, 0, 0,0,0,0,0), 

2 . Gauge field mode (R): <5 = —s/s and v = (0,0, 0, 0,1, 0, 0), 

3. Constant shift of the leading helix parameter wq (M): 5 = — ^/s and 

V = (0,0, 0,0, 0,1,0), 

4. Charged scalar mode (if real: 6 - (R), <5+ (I) if rUp > 0 or (R) if < 0, (M) if 

"ip = 0): = -| (1 T Y^l + and v = (0,0,0,0,0,0,1). 

5. Combined gauge and helix field mode: = g ^—9 ± Vl + 120 k^^ and 

v± = ^0, 0, 0, 0, ’ "''^hh (5_ < 0 always (R), and <5+ < 0 (R) (5+ > 0 

(I)) for |k| < « 0.817 (|k| > ~ 0.817), 

6. Combined matter and geometry mode: = g (—5 ± \/l45) (5+ (I), 6 - (R)) and 

v+ = (l4 - V 145 ) , -1, -1 (3 + V 145 ) , 1, 

10 (3 (116105 + 9643\/l45) + I081\/l45 + 13103) 

17 (15 (26309 + 2185\/l45) _ 70009\/l45 - 843005) ’ 

4 (-12 (7720 + 641\/l45) + 8071\/l45 + 97211) 

17 (3 (4105 + 341\/l45) _ 2185\/l45 - 26309) ’ 

(l4 + v/l45) , -1, -1 (3 - v/l45) , 1, 
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10 ((348315 - 28929/145) / - 1081/145 + 13103) 

17 (15 (2185/145 - 26309) k? - 70009/145 + 843005) ’ 

4 ((92640 - 7692/145) / + 807l/l45 - 97211) \ 

17 (3 (341/L45 - 4105) k? - 2185/L45 + 26309) ’ ) ’ 

7. Combined geometry mode: / = g (—5 ± /l85) (<5+ (I), 5- (R)) and 

v± = (71 T 5/1^) , (31 T 2 / 1 ^) , (31 T 2/1^) , 1 , 

^ (205 T 17 / 1 ^), ( 31 T 2 / 1 ^), 0 ^ , 

8. Combined matter and geometry mode (M): <5 = 0 and v = (0, —1, 2,1, 2,4,0), 

9. Combined matter and geometry mode (R): <5 = —1 and v = (6, —1, 2,1, 5,4,0). 

Here modes 1,6,7,9 can contribnte to the temperatnre mode. Besides the obvions pairs 
(4,5,6,7) which snm np to —/a (4,6,7) and —3 (5), there seem to be single modes (3,8,9) as 
well. 

C.3 Condensed Fixed Point 

For convenience, we switch back to the Ansatz for the pertnrbations (4.3). Within the 
system (2.2) inclnding the charged scalar, onr insnlating geometry with charged scalar hair 
(4.1) has the following radial pertnrbations: 

1. Combined blackening factor, gange field and charged scalar mode (R): 5 = —/s and 

V = (—1, 0, 0, 0, —1, 0, 5) (corresponds together with mode 2 below to mode 1 and 2 
in Appendix C.2), 

2. Combined blackening factor, gange field and helix mode (R): <5 = —/a and v = 

, 0,0,0, (corresponds together with mode 1 above to mode 

1 and 2 in Appendix C.2, pQ ^ 0 not obvions), 

3. Constant shift of the leading helix parameter tco (M): <5 = —/s and 

V = (0,0, 0,0, 0,1,0), 

4. Charged scalar modes: 

Constant shift in po (M): <5 = —/s and v = (0, 0, 0, 0, 0,0,1) (corresponds to 5+ of 
mode 4 in Appendix C.2), 

IR relevant scalar mode: <5 = —3 and v = (0,0,0,0,0,0,1) (corresponds to of mode 
4 in Appendix C.2 for = 0, as for this mode Sp ~ 

5. Combined gange, helix and charged scalar mode:^® 

= g ^—9 ± /l + 120 k^ + 20//9g^ and 

^®Compared to [26], there is a shift of 2/3 in this exponent because the perturbations of the matter fields 
are written differently with respect to the background. 
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V± = 


n n n n ('?^/’0-4)(-1 ±a/1+120k2+20(?2p 2) 3(i}2p§-4)(-9±y^l+120)?+20g^) 

48k 2 10K2(_4+6K2+52p2) 


with 


(5_ < 0 always (R), and (5+ < 0 (R) (5+ > 0 (I)) for 1 + 120 k^ + 20 g^/9Q < 81 

(1 + 120^2 + 20gVo >81), 

6 . Combined matter and geometry mode: = g (—5 ± \/l45) (<5+ (I), 6 - (R)) and 


v+ = 



yi45 1 
6 ’ 6 




i(l4 + yi45), 


30 (71 + 7145) 7 + 5 (71 + 7145) q^pl + 1787145 - 2050 
6 (6 (77145 - 115) ^2 + (77145 - 115) q^pl - lOVuE + 514) ’ 

2 (67 ((1295 + 7145) - 4 (1565 + 9l7T45))) 

3 (7p2 _ 4) (0 (477145 _ 335) ^2 (477I45 _ 335) 7^2 _ 2427145 - 46) 

2 (7p2 _ 4) 71295 + 7145) q^pl - 2 (5209 + 2637145)) 

3 (g2p2 _ 4) (0 (4771^ _ 335) ^2 (4771^ _ 335) q2p2 _ 2427145 - 46) ’ 

2 (6 (417145 - 1985) 7 + (417145 - 1985) q^pl - 6147145 + 10262) \ 

18 (205 + 237145) k2 + (015 + 697145) q^pl - 7987145 - 8922 J 

30 (7l45 - 71) 7 + 5 (7l45 - 7l) q'^pl + 1787145 + 2050 
6 (6 (115 + 77I45) k 2 + (115 + 77145) g2p2 _ 2 (257 + 57145)) ’ 

2 (6k2 ((7145 - 1295) - 3647145 + 6260)) 

3 {q'^pl - 4) (6 (335 + 477145) ^2 + (335 + 477I45) q'^pl - 2427145 + 46) 

2 ((7/j2 - 4) ((7145 - 1295) - 5267145 + 10418)) 

3 {q'^pl - 4) (6 (335 + 4771^) ^2 + (335 + 477!^) gVo “ 2427145 + 46) ’ 
2 (6 (1985 + 4171^) 7 + (1985 + 4l7l^) q^pl - 2 (5131 + 3077!^)) \ 

18 (2371^ - 205) k2 + (697145 - 615) - 79871^ + 8922 j ’ 


7. Combined geometry mode: = g (—5 ± 7l85) ((5+ (I), 6 - (R)) and 

''-^= (‘7 7 "-^) 4 7 +^) 7 

10 (3 (35 + 37185) + (67 + 57185) 7po “ 22 (l3 + 7l85)) 

30 (21 + 7185) 7 + 5 (21 + 7185) 7Po “ 4 (355 + 237185) ’ 

(<77o - 4) (5 (87 + 77185) q^pl - 5247185 - 6940) - 120^2 (Sg^pg + 7^135 + 37 ) 

(g 2 p 2 _ 4 ) (30 (31 + 37 M) ^2 5 (31 3 ^ 1 ^) ^ 2^2 _ 4 (005 53 ^ 1 ^)) 

2 (60 (83 + 97185) 7 + 5 (1057 + 8l7l85) g^Po “ 176 (85 + 77185)) \ 

30 (479 + 357185) ^2 + 5 (479 + 357185) g 2 p 2 _ 4 (9225 + 6777185) J 

= ( 7 ('^- 777 -^) 7 ('^'^- 7 - 

10 (3 (37185 - 35) 7 + (57185 - 67) q^pl - 22 (7l85 - 13)) 

30 (7l85 - 21) 7 + 5 (7l85 - 2l) g2p2 _ 9271^ + 1420 ’ 
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(g2p2 _ 4) (5 (77185 - 87) q^pl - 52471^ + 6940) - 120k2 - 87) 

(gVo - 4) (30 (37185 - 3l) + b (371^ - 3l) q'^pl - 21271^ + 2660) 

2 (60 (971^ - 83) k 2 + 5 (8l7l^ - 1057) q'^pl - 176 (771^ - 85)) 

30 (357M - 479) /c2 + 5 (3571^ - 479) q'^pl - 270871^ + 36900 

8 . Combined matter and geometry mode (M): <5 = 0 and v = (0, —1, 2,1, 2,4,4), 

9. Combined matter and geometry mode (R): 5 = —1 and v = (6, —1, 2,1, 5,4,4). 

Here modes 1,2, 6, 7)9 can contribute to the temperature mode. Besides the obvious pairs 
(4,5,6,7) which sum up to —s/s (4 after taking into account the different Ansatz for the 
fluctuations, 6, 7), and —3 (5), there seem to be single modes (3,8,9) as well. 



D Numerical Method for Background and Fluctuations 


From a numerical perspective, we have to solve a boundary value problem for a set of 
coupled, non-linear, ordinary differential equations. This is done using a shooting method 
consisting of the following steps: 

1. Choosing an initial guess of horizon parameters, the asymptotic horizon expansion 
is used to set up initial conditions rp + 5 with 5 numerically small. The horizon 
radius rp can be set to one by a radial rescaling {c.f. the scaling symmetries discussed 
below). 

2. The equations of motion are integrated numerically between rp + S and ^ rp using 
Mathematica’s numerical integrator NDSolve. 

3. The difference between the numerical solution and the desired boundary values at 
is read out. 

4. The integration between the horizon and the boundary is iterated. Using Broyden’s 
method [71] as a root finding algorithm, the horizon parameters for which the numer¬ 
ical solutions satisfies the boundary conditions are determined. 

5. The boundary parameters are determined by matching the asymptotic boundary ex¬ 
pansion to the numerical solution. 


Alternatively, some of the boundary conditions can be imposed making use of the following 
scaling symmetries of the equations of motion: 


(I) 

r —>■ yr. 

t —^ 

(11) 

X —)• 7X, 

p —)• 

(III) 

(y,z) ^{y,z), 

w —)■ 


U 

e 2 ^i 

g2l)2,3 




e 2 i'i 
9 ’ 

T 

q2v2,3 


T 


a —)■ 70, 


(D.l) 
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Using scaling (II) and (III), we can set e^*'^ = for r ^ and either e^’^^ = or 
q 2 v 3 — j .2 ^ay, the scaling symmetries reduce the number of boundary conditions 

which have to be imposed on the numerical solution by means of the shooting method. They 
are in particular useful for finding a first solution to the equations of motion. However, the 
scalings (II) and (III) change the value of p and the asymptotic value of the helix field w. 
Once a branch of solutions is found, it is therefore more convenient to use the shooting 
method as explained above to satisfy all boundary conditions. Once the backgrounds are 
generated, we use the following numerical procedure to calculate the conductivity for a 
given solution: 

1. Using Mathematica’s NDSolve, three linearly independent sets of solutions are con¬ 
structed by numerical integration between + S and (with <5 <C 1 and ri, ^ Vh). 

2. The three sets of solutions are linearly combined to a solution for which Q and /i23 
have a vanishing source and which satisfies the (arbitrary) normalization condition 
A[n) = 1. 

3. The boundary expansion modes of the fluctuations are determined by matching the 

asymptotic expansion valid for r to the numerical solution. 

Given that the work of [26] was not very explicit about the number of free UV and IR 
parameters in the zero temperature RG flows of our model, let us conclude this section by 
a more detailed explanation of the situation: As noted above, due to the broken conformal 
symmetry (2.37) and the five asymptotic conditions 

U(r)~r^, •Uj(r)~lnr, i = l,2,3, Jp = 0, (D.2) 

we expect that the zero temperature RG flows are labeled by the two parameters (A, p). In 
the IR, the conditions (D.2) are fixed by the five IR parameters (po) '*^20 and '*^30; and the two 
coefficients of the two IR irrelevant modes (4.4)). We are hence left with a two-parameter 
family of solutions, labeled by the two parameters (A, p). The vacuum expectation values 
in the UV are then fixed by the requirement of regularity in the IR, i.e. by the vanishing of 
the coefficients of the IR relevant perturbations around the condensed geometry (the modes 
under point 1,2,4,5,6,7,9). Dropping the charged scalar both in the UV and in the IR, the 
same argument applies to the uncondensed insulating solutions. The metallic AdS2 x 
solutions however are unique {i.e. a zero-parameter family) due to the vanishing of the helix, 
A = 0, as well as the charged scalar, and hence the restored conformal symmetry (2.37). 
This is matched in the IR by two free parameters t>io, V 20 (^20 = '^^3o) the coefficients 
of four IR irrelevant perturbations (modes 7,8,9,10 in Appendix G.l), and hence there 
are no free parameters left (except of p) in this holographic RG flow. In particular the 
vanishing of mode 8 and 9 in Appendix G.l sets the charged scalar and the helix field to 
zero, respectively. An important outcome of this whole discussion is that the helix pitch p 
should not be counted as an independent UV parameter, as it is not a source in the usual 
quantum field theoretic sense, but rather a parameter of the boundary geometry. 
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E Holographic Renormalization and Operator Mixing 


In order to calculate the retarded Green function, the on-shell action of the fluctuations is 
evaluated. As for the background, the on-shell action can be reduced to a boundary term 
by partial integration since the bulk term vanishes upon use of the equations of motion. 
Fourier transformation of the fluctuation fields, 

A{t,r) = I ^e--M(a;,r) (E.l) 


and similarly for the remaining fluctuation fields, results in an action bilinear in the Fourier 
modes. It has the form 


*5, 


dwdrd^x 




(E.2) 


where $ is the vector of fluctuation fields and Mj are matrices of r-dependent functions 
containing the background fields. The derivatives of modes that are evaluated at (—w) can 
be eliminated by partial integration. Upon use of the equations of motion the bulk term 
vanishes and Sq reduces to a boundary term. The horizon contributions to this boundary 
term are discarded following the prescription of [47]. To regularize divergences, we introduce 
an ultraviolet cutoff r^, which will be removed eventually. The full on-shell action of the 
fluctuations Sos is given by the sum of Sq, the Gibbons-Hawking term 

5gh = 2 y' dt d^x (E.3) 


expanded to second order in the fluctuations, and appropriate counterterms ensuring 
that Sos = Sq + S'gh + 'S'ct is finite in the limit rf, —)• oo. The counterterms needed to make 
Sos finite are 


si = / dtd^x^/^- 


1 


1 , 


1 


-6 - -R^ + log(rfe) ( -FabF^^ + -W'^'^Wab - ^R^^’R 




■ (E.4) 


Here 'jab denotes the induced metric at r = rb, R-y is the Ricci scalar of the induced metric, 
and R'y^ab is the Ricci tensor. The studies [48, 72] were taken as references in finding possible 
counterterms. In order to extract the retarded Green function, we can either switch off the 
sources of the physical fields G and F by hand, or use the holographic operator mixing 
method [48, 73] as explained in the following. Arranging the physical fields in a vector, 
$(phys) _ the terms of the on-shell action containing these fields can be written 


as 






J r=rb 


(E.5) 


Here Ma and Mb are the two matrices 


U 


/ qV2+V3-Vi 

0 

0 1 


/- 

e-'"i+’'2+”3 u]2 log(r) 

K.W 



2 VU 

4 

0 

e»'i+'“2-"3 

0 

and 


KW 


. 

U}^ 


4 

* 


V 0 

0 

gl)l-|-3D2-t'3 


1 

0 

* 
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respectively. The entries marked with an asterisk correspond to rather longish functions 
of r containing the background fields. Their precise form is not needed since we are only 
interested in the retarded Green function corresponding to A. Numerically, we can con¬ 
struct three sets of linearly independent solutions and by integration of 

the equations of motion with three linearly independent initial conditions. Arranging the 
solutions in a matrix, 


H = 


( I I 

^phys ^phys 

V I I 



(E.6) 


the matrix of the retarded Green function corresponding to {A, Q, J-) can be calculated as 


Gr{oj) = 2 n)H'{uj, n)H n) - Mb{uj, rb)) 




(E.7) 


A derivation of this result can be found in [72, 73]. The basic idea consists in constructing 
sets of solutions to the equations of motion with each set sourcing only one particular 
fluctuation on the boundary. Finally, a generalization of the prescription of [47] to the case 
of multiple fluctuation fields results in the above formula for the retarded Green function. 
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